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INTRODUCTION 


Landing delays are currently a problem at major air carrier airports, and most 
forecasters agree that they will get worse by the end of the century. It is 
anticipated that some types of delays can be reduced by an efficient high 
speed runway exit system, allowing the increased approach volumes necessary to 
assist congested airports. 

This report defines a computerized Probabilistic Runway Turnoff Model 
developed in Ref. 1, which analyzes exit utilization and operational 
characteristics for varying aircraft types. 

To achieve an increase in arrivals at congested airports, methods of 
minimizing both runway occupancy time and its related standard deviation must 
be achieved. Reliability is the major factor — thus it was established that 
one miss for every 10,000 turnoffs per week per airfield (2 runways) or 3.75 
standard deviations above average, is an appropriate goal. When considering 
congested existing airports only experiencing approximately 2,000 weekly 
landings per runway, compared to the goal of 10,000, the 1/10,000 reliability 
looks extremely attractive. (See Table 1) 

With this level of reliability, separation distances between aircraft could be 
decreased to provide an increase in density of arrivals. 
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Table 1 

RELIABILITY OF RUNWAY OPERATIONS 


Number of 
Weekly Landings 
Per Runway 

Probability 
of One Missed 
Turnoff Per Week 
Per Runway 

No. of 
Std Devs. 
From the Mean 

Example Airport 
Experiencing This 
Activity per Runway 

5,000 

.0002 

3.75 

N/A 

4,000 

.00025 

3.45 

N/A 

1,700 

.00059 

3.23 

LaGuardia 

1,200 

.00083 

3.15 

LAX 




Denver 


Ref. 2 


Total aircraft movements (52 wks/yr x 2 (approaches/departures) x number 
of runways used for departures.) 
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PROBABILISTIC RUNWAY TURNOFF MODEL 


PROGRAM FUNCTION 


The computer probabilistic model is comprised of two parts; 1) the time 
required from threshold to start of exit with a probability determination of 
an exit velocity and 2) a subroutine of time required in the turnoff to clear 
the runway using an optimized path. The times determined by each part are 
added together to yield the total runway occupancy time which, 
probabilistically, will be a unique value less than 40 seconds. The flow 
chart in Table 3 outlines the inputs, the major steps, and the iteration loops 
of the probabilistic runway turnoff model. 

RUNWAY OCCUPANCY TIME 


Runway occupancy time is defined as the average time plus the number of 
standard deviations for a required probability. Therefore, a maximum 
occupancy time limit, not an average, of 40 seconds or 90 landings per hour, 
was incorporated with the following rational. 

If an average time was used in the analysis, assumed to be normally 
distributed, there would be a substantial number of aircraft with occupancy 
times approaching 30 seconds. However, the aircraft must be spaced on final 
approach with the intent of a 40 second occupancy time. Arrival separations 
and runway occupancy times are mutually exclusive events. The Air Traffic 
Controller is not capable of determining the amount of time any given aircraft 
will spend on the runway. Therefore, it was determined that the only real 
alternative with which to evaluate runway exit design would be to use a 40 
second maximum. This time measurement is from the threshold to the point the 
aircraft has completely cleared the runway surface during exiting. 

TOUCHDOWN SPEED AND LOCATION 


Touchdown speed and location are assumed normally distributed even though 
existing data indicate the log-normal distribution is a better fit of 
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touchdown location. The time from touchdown until start of deceleration is 
assumed normally distributed, and the speed is assumed constant from crossing 
the threshold to the start of deceleration. 

BRAKING DECELERATION 


Since most aircraft are capable of decelerating faster than is normally 
practical, a conservative deceleration rate is used. 

The deceleration process has a number of factors contributing to the total 
deceleration force. As Figure 1 shows, the major force is from the 
application of brakes, followed by aerodynamic drag at higher speeds and then 
reverse thrust. However, during deceleration, modern aircraft hold a constant 
deceleration rate controlled by avionics equipment which controls brakes at 
increased levels after reverse thrust is shut down and aerodynamic drag has 
decreased. This auto-braking system allows different combinations of braking 
versus thrust reversal to be utilized to account for different weather 
conditions. Therefore, a constant deceleration rate is maintained down the 
runway to the beginning of the exit path. Because all of the pilot’s 
attention is required during his exiting procedures, reverse thrust is shut 
down before exiting rather than during exiting and deceleration is decreased 
through the path, helping maneuverability performance. 

This deceleration rate named XA, Table 2, is an input to the model which 
allows for changes of deceleration rates to the user’s choice. 

PROBABILITY EXIT SPEED 

Speeds and distances from the exit entrance are assumed normally distributed. 

A 50-percentile cumulative normal distribution speed, named YA, initializes 
the program. A higher YA speed (YA + YS) yields a high probability of 
combined variables achieving 40-second maximum runway occupancy time, while a 
lower YA speed (YA - YS) gives a lower probability. The approach to finding 
the 99 percent to 100 percent probability speed of exiting with a 40-second 
maximum runway occupancy time is achieved. 
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If the 40-second criterion is not met, the time is proportioned based on 40 
seconds, and repetitive substitutions of deceleration and exit to threshold 
distances (CDIST) are made until the 40 seconds maximum is met. 


CF6-6D 


o 

o 

o 


LLf 

o 

DC 

o 



FIGURE 1. SUMMARY OF AIRPLANE FORCES DURING LANDING GROUND RUN 
MODEL DC-10 SERIES 10 
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Table 2 


PROGRAM LABELS 


Parameter 

Average 

Standard 

Deviation 

For aircraft traveling 
one standard deviation 
below the average speed 

Touchdown location 
(distance from threshold) 

UA 

us 




Touchdown speed 

VA 

VS 




Time from touchdown until 
start of deceleration 

WA 

ws 




Deceleration rate 

XA 

xs 




Probability speed 

YA 

YS 




Speed at which reverse 
thrust stops 

ZA 

SDB 


SB* 


Speed at exit of study 

ASC 

SDC 


SC' 


Speed at point where 
deceleration starts 
For aircraft traveling 

ASA 

SDA 


SA' 




Average 

1 

For aircraft traveling 
one standard deviation 
below the average speed 

Distance from threshold 
where deceleration starts 


A 


A' 


Distance from threshold 
where reverse thrust stops 


B 


B' 


Distance to threshold of 
exit being studied 


C 


C' 


Average runway occupancy 
time at A 


TA 


TA' 


Average runway occupancy 
time at B 


TB 


TB’ 


Average runway occupancy 
time at C 


TC 


TC' 
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Distances UA, A, B, and C in Figure 2 are used through the flow chart, Table 3, 

and the reader must be familiar with them to adequately understand the model 
logic. 



FIGURE 2. RUNWAY OPERATIONS 


ALTERNATE PROGRAM MODES 


When initially developed, the model was programmed to solve the 40 second 
occupancy time criterion via the Mode 1 iterative process. Since iteration 
leaves too many variables, accuracy in every case becomes difficult. 

Therefore, Modes 2 and 3 were introduced to allow the user to choose his own 
values for YA. 

Mode 2 allows the user to initiate a YA, then the program reduces YA by a 
value of 1.0, for 15 iterations through the probabilistic model and prints the 
values in descending order. This becomes useful when evaluating the impact of 
YA on exiting probability and runway occupancy time for a closed data set. 


Mode 3 is the quickest method of obtaining direct solutions for a given YA 
value. After each completed run, the user is asked to submit a new YA value 
if not satisfied with the probability. When satisfaction is reached, the user 
enters 999 to print a completed list of the previously entered YA values and 
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TABLE 3 

FLOWCHART - PROBABILISTIC RUNWAY 
TURNOFF MODEL 
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TABLE 3 (cont.) 








TABLE 3 (coni.) 


EXPATH SUBROUTINE 
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their probability values. If the probability is less than 1.000 the YA value 
must be increased. However, if the probability is greater than 1.000, YA must 
be decreased to a probability of less than 1.000 to determine the breaking 
point of 1.000 and 0.999. This breaking point identifies the specific data 
set to be used in the input data. Tables 4-6 exemplify the three types of 
modes described above. 


TABLE 4 

MODE 1 EXAMPLE 


ENTER: 5 FDR TERM! URL PRINTOUT 
6 FOR HR Pit CORY PRINTOUT 
5 

Rs 1 FDR ST HUP RED PRINT 

2 FOR INTERMEDIATE VALUES 

3 FOR FINAL VALUES ONLY 
^1 » 

R: 1 TO SOLVE FOR 40 SEC OCCUPANCY 

2 TO SOLVE FOR MAXIMUM PROBABILITY 

3 FDR USER ENTRY OF YA VALUE 

1 

«" ER U,ITIRL Vfi SPEED REF. TERPS CATEGORY E AIRCRAFT 

O i.l M 

CALCULATED RESULTS 


CRITICAL VALUES: 

AS CSV H> 

276.22 


PTIME 


71 


26-4 . 0 1 


6. 06 


CD I ST 
34 0 0 . 0 0 


XDIST 


1634.71 


PROB 

AV6 TIME 

STD DEV 

1 - 0 00 0 

£4 . 76 17 

£.3551 

0. 0000 

c’4. 7617 

£.3551 

0. 0000 

£4.7617 

£.3551 

0. 0000 

£4 • 76 1 7 

£.3551 

0. 0000 

£4. 7617 

£ . 355 1 

0. 0000 

£4. 7617 

2.3551 

0. 0000 

£ 4 . 761 7 

£ . 355 1 

0. 0000 

£ 4 . 7617 

£.3551 

0. 0000 

£4 . 7617 

£.3551 

0. 0000 

£4. 7617 

£.3551 

0. 0000 

£ 4 . 7617 

£.3551 

0. 0000 

c.’ 4 • i' 6 1 i*’ 

£.3551 

0. 0000 

£4 . 76 1 7 

£.3551 

0. 0000 

£4.7617 

£.3551 

0 . 0 0 0 0 

£4.7617 

£.3551 


— STOP 


XT I ME 

TTIME 

Vfi 

33. 9466 

•4 0. 01 04 

3 0 0 . 0 0 0 

33. 9466 

4 0 . 0 1 04 

£54. 5996++ 

33. 9466 

4 0. 01 04 

£16. 0699 + 

33. 9466 

4 0. 01 04 

133.3711+ 

33. 9466 

4 0 . 0 1 04 

1 55 . 6£ 07* 

33. 9466 

4 0. 01 04 

1 3 £ • i j 6 9 9 * 

33. 9466 

4 0. 01 04 

1 1£. 0S3£+ 

33. 9466 

4 0. 01 04 

95. 1 £ 1 1 ♦ 

33'. 9466 

4 0. 01 04 

SO*. 7£60+ 

33. 9466 

4 0 . 0 1 04 

63. 5*094 + 

33. 9466 

4 0. 01 04 

53. 1415+ 

33. 9466 

4 0 . 0 1 04 

49. 34 £7+ 

33. 9466 

4 0. 01 04 

41.3754+ 

33. 9466 

4 0. 01 04 

35. 533£+ 

33. 9466 

4 0. 01 04 

30. 1601+ 
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TABLE 5 


MODE 2 EXAMPLE 


ENTERS 5 FDR TERMIHRL PRINTOUT 
6 FDR HARDCOPY FRIHTOUT 
5 

R: 1 FDR STANDARD PRINT 

£ FDR INTERMEDIATE VALUES 
3 FDR FINAL VALUES DULY 

J 

Rs 1 TD SOLVE FDR 40 SEC OCCUPANCY 
£ TD SOLVE FDR MAXIMUM PROBABILITY 
3 FOR USER ENTRY DF YA VALUE 

£ 

ENTER INITIAL YA SPEED 
300 

TERRS CATEGORY E - LOCKHEED F-104 


NASA HIGH SPEED EXIT ANALYSIS 
INPUT SUMMARY 


PROBABILISTIC MODEL 


TOUCHDOWN LOCATION 
TOUCHDOWN SPEED 
TIME-TOUCHDOWN TD DECEL 
DECELERATION RATE . 

5 O’: PROBABILITY SPEED 
BIST: THRESHOLD-START DF TURN 
SPEED AT MAJOR DECEL REDUCE 
NDSE-WINGTIP DIST 


AVERAGE 

1500. 

370. 

5. 

6 . 0 0 
3 0 0 . 

84 0 0 . 

1 5 0 - 
35. 


STD DEV 
3 0 • 

£5. 

1. 

£ . 0 0 

1. 


UNITS 

FEET 

FT/SEC 

SECONDS 

FT/SEC++2 

FT/SEC 

FEET 

FT/SEC 

FEET 


AIRCRAFT TRACKING PROGRAM 


A I RCRAFT CHARACTER I ST I CS 

WEIGHT <LBS> YAW INERTIA WHEELBASE LOAD OH MAIN GEAR 

££000. 0 137776. 0 15.0 9£. 3 


FLIGHT CONDITIONAL PARAMETERS 
EXIT VELOCITY MAX EXIT ANGLE 
£76. £25 £0. 000 


MAX Y MU DECEL RATE 

0. £00 — £. 000 


SYSTEM WIDTHS 
RUNWAY 
150 . 00 


WINGSPAN 

£1.90 


CLEAR DISTANCE 
85. 95 


SCRUB CHARACTER I ST I CS 
Y MU RADIUS 


0. 0000 

0. 0200 


0. 1 E+32 
50. 0 
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TABLE 5 (contd.) 


CRLCULfiTEIi RESULTS 


CRITICAL VALUES: 


ns csv 

fiSC 

FT I ME 

CD I ST 

XD I ST 

d \ t> • d d 

d 6 4 • IJ 1 

£>. 06 

84 00. 0 0 

1634 


PRDB 

fiV6 time 

STD DEV 

XT I ME 

TTIME 


Yfl 

1 . 0 0 0 0 

34 . 76 1 7 

£.3551 

33. 9466 

4 0. 01 04 


00. 00 C 

1 . 0 0 0 0 

84.7617 

£.3551 

33.9466 

4 0 . 0 1 04 

339. 

0 0 0 »> ♦ 

1 . 0 0 0 0 

84. 7617 

£.3551 

33.9466 

40. 01 04 

898 . 

0 0 0 0+ ♦ 

1 - u 0 0 0 

d 4 • i*’ t:* 1 »* 

£.3551 

33 . 9466 

4 0. 01 04 

897. 

0 0 0 0+ * 

1 . 0000 

d 4 • t 6 1 * 

£.3551 

33.9466 

4 0. 01 04 

896 • 

ij 0 tf 0 ♦ ♦ 

1 . 0 0 0 0 

84.7617 

£ . 355 1 

3 3 • 946 6 

4 0 . 0 1 04 

895. 

0 0 0 0 ♦ ♦ 

1 . 0 0 0 0 

84. 7617 

£.3551 

33.9466 

40. 01 04 

894. 

00 0 U+-* 

1 . 0 0 0 0 

d 4 • i 6 1 * 

-:*cc i 

33. 9466 

4 0. 01 04 

893. 

00 0 0** 

1 . oooo 

84. 7817 

8.3551 

33. 9466 

4 0. 01 04 

898. 

0 0 0 0* v 

fi. 9999 

84. 7617 

8.3551 

33. 9466 

4 0. 01 04 

891 . 

i j i j i j 0 ♦ ♦ 

O. 9975 

84. 7986 

3.3551 

33.9335 

4 0. 0473 

890. 

0 00 1> - 

ft.9734 

84.9313 

8 • 355 1 

34. 1168 

4 0. 1799 

839 . 

000 

0 • 87 6a 

85. 3464 

8.3551 

34.5313 

4 0.5951 

838 . 

Ij 0 Ij Ij -r * 

0.7880 

85.9039 

8.3551 

35. 0938 

41.1 576 

887 

. ijijij0*< 

0.6196 

d 6 • ubc'4 

8.3551 

35.4473 

41.5110 

886. 

Ij Ij Ij Ij* ♦ 


— STDP 
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table: 6 

MODE 3 EXAMPLE 


ENTER: 5 FOR TERMINAL PRINTOUT 
6 FDR HARDCOPY PRINTOUT 

5 

1 FDR STANDARD PRINT 

2 FOR INTERMEDIATE VALUES 
Si FOF! FINAL VALUES ONLY 

R: 1 TO SOLVE FDR 40 SEC OCCUPANCY 

2 TO SOLVE FOP MAXIMUM PRODAB1L ITY 

3 FOR USER ENTRY OF YA VALUE 


ENTER INITIAL YA S 
6000 

PROB = 1.0000 

ENTER NEW YA SPEED? 

F'l-IED 

TO 

END 

£0 

I; = 0 . 0 0 0 0 

ENTER NEW Yfl 

speed; 

999 

TO 

END 

•400 

DB = 1 . 0 0 0 0 

ENTER NEW Yfl 

speed; 

999 

TO 

END 

£ 0 0 

B * « 0 . 0 0 0 0 

ENTER NEW Yfl 

speed; 

9 9 9 

TO 

END 

300 

B = 1.00 0 0 

ENTER NEW Yfl 

speed? 

9 9 9 

TO 

END 

270 

B = 0.0029 

R1TER NEW YA 

speed; 

999 

TO 

END 

295 

PROB = 1.0000 

ENTER NEW YA 

speed; 

999 

TO 

END 

290 

ENTER NEW* YA 

speed; 

9 9 9 

TO 

END 

294 

B = 1 . 0 0 0 0 

ENTER NEW YA 

speed; 

9 9 9 

TO 

END 

E = 1 . 0 0 0 0 

ENTER NEW Yfl 

speed; 


TO 

END 

292 

B = 1 . 0 0 0 0 

ENTER NEW Yfl 

SPEED? 

9 9 9 

TO 

END 

291 

E = 0.9999 

ENTER NEW Yfl 

speed; 

999 

TO 

END 


99 ? 


REF. TERPS CATEGORY E AIRCRAFT 
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TABLE 6 (contd.) 

CALCULATED RESULTS 


CRITICAL VALUES: 

ASCSV A SC PTIME CD I ST XDIST 


d l' - cl 

cl c‘64 . 

01 

6. 06 

8400. 00 

1634.71 

PRDB 

fiVG TIME 

STD DEV 

XT I HE 

TTIME 

V ft 

1 . 0000 

£4 . 7617 

£.3551 

33.9466 

40.01 04 

6000. 0000++ 

0. 0000 

£4 • 761 7 

£.3551 

33. 9466 

4 0. 0104 

£0. 000 0^ 

1 . 0000 

£4 « 76 1 7 

£. 3551 

33. 946*6 

40. 01 04 

400. 0000** 

0. 0000 

£4 • 76 1 7 

£.3551 

33. 9466 

4 0. 01 04 

£00. 0000+ 

1 . 0000 

£4 • 76 1 7 

£.3551 

33. 9466* 

4 0. 0104 

300. 00 0 0*+ 

0. 00 £9 

£4.80£1 

£.3551 

33. 9869 

40. 0507 

£70. 0000++ 

1 . 0000 

£4.7617 

£.3551 

33.9466 

40. 0104 

£95. 0 0 0 0 ♦ ♦ 

0. 9975 

£4.7986 

£.3551 

33.9335 

4 0. 0473 

£90. 00 00** 

1 . 0000 

£4 . 76 1 7 

£.3551 

33. 9466 

4 0. 01 04, 

£94. 00 00+* 

1 . 0 0 0 0 

£4.7617 

£.3551 

33. 9466 

4 0. 01 04 

£93. 00 0 

1 . 0 0 0 0 

£ 4 . 7 c. 17 

£.3551 

33.9466 

4 0. 01 04 

£9£. 0 0 0 0* ♦ 

0. 9999 

£4.7617 

£.3551 

33. 9466 

4 0. 0104 

£91 . 000D+* 


— ST DP 


EXIT PATH (TRACKING) SUBROUTINE 

The exit path subroutine is essentially the same program as the one developed 
in Ref. 1 with the following exceptions: 

Deceleration during the turn has been added as an input. 

- Wing tip location has been added to tell when the aircraft has 
cleared the runway. 

Passenger comfort limitations (G and J_ a „) have been 

max max 

internally set to inifinity so that they are not limiting. 

SUBROUTINE DEVELOPMENT 


The exit path optimization subroutine flow chart is shown in Table 3. 
Basically, the subroutine computes the nose gear / u consumed In overcoming 
scrubbing and centrifugal forces, applies the net [l remaining to the yaw 
inertia to find the rate at which the turning radius is shrinking, and then 
integrates this rate of turn radius change along with the other motion terms. 
The integration step size is 0.001 seconds. Data output at print intervals is 
just prior to the integration equations. Figure 3 shows input nose gear ys, 
required to overcome scrubbing forces. 
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I I I I I I 

0 100 200 300 400 500 

(METERS) 

TURNING RADIUS 


FIGURE 3. NOSE TIRE SIDE M's VERSUS TURNING RADIUS DURING SLOW-SPEED TURNS 


LANDING GEAR FRICTION CHARACTERISTICS 


The ability of the nose gear to control the direction of the aircraft is a 
direct function of (1) its vertical load, (2) the side load it must develop to 
accomplish a particular maneuver, and (3) the ground coefficient of friction. 
For simplicity, the demand placed on the nose gear is expressed as the ratio 
where : 


Cnq = Side Load on Nose Gear 

Vertical Load on Nose Gear 
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This allows direct comparison to the available ground coefficients of friction 
which are : 


Dry pavement: 
Wet pavement: 
Packed snow: 
Ice: 


0.8 

0.4 (varies with speed) 
0.2 
0.1 


The demand on the nose gear comes from three sources : 
1. Centrifugal force: 


C ay » V 2 /(g c R) 

where 

V « velocity 

6 C 55 acceleration of gravity 

R - instantaneous turn radius 


2. Rotation inertia resistance: 


C I = Ioc/(g c ) (W) (L w ) (XM/100) (1-%M/100)) 


where I 
oc 


rotational inertia of aircraft 

rotational acceleration in radians/sec 2 = -(V) (R')/R^ 


where R' 
W 
Lw 


is the rate of change in turn radius 
aircraft weight 
aircraft wheelbase 

percent of gross weight supported by the main landing gears 


3. Scrubbing resistance: 


Figure 3 shows pi's for five airplanes. The aircraft's resistance to 
being turned is dependent upon the type of gears used and their location. 
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Duals are very easy to turn, whereas a configuration of two wing-mounted 
gears and two body-mounted gears 3.0 m aft is the most difficult of the 
commercial jets to turn. 

Total maneuvering demand on the nose gear then is : 

C = Ca„ +C T + C c = JU's = LL 
n S y I s s “ r 1 scrub 

^ ^ n g 0.2 is considered to be a practical limit for maneuvering at speeds 
of 59-89 ft/s. This would leave the pilot with a margin of 4 against skidding 
on dry surfaces and a margin of 2 against skidding on wet surfaces. Pilots 
would probably refuse to consistently use margins less then 2. 

Note that low weight, aft center of gravity conditions are the most critical. 
Aft center of gravity results in the least weight to the nose gear. Both 
rotational inertia and main gear turning resistance remain fairly high at low 
weights. The listing in Table 7 is typical of the airplane data required by 
the program. 


TABLE 7 

AIRPLANE DATA FOR RUNWAY TURNOFF PROGRAM 


Aircraft 

Weight 

LB 

Yaw Inertia 
10 6 Lb In2 

Wheelbase , LW 
FT 

Wt Supported 
by Main Gear 
% 

DC-9- 10 

65,372 

4546 

43.6 

95.4 

DC-9-80 

102,418 

17827 

72.6 

96.8 

DC-8-63 

173,630 

51545 

77.4 

96.7 

DC-10-10 

290,850 

75813 

72.6 

94.2 

747-100 

451,277 

194627 

84.0 

96.4 


18 











PASSENGER COMFORT 


Equations and logic are included in the subroutine for additionally limiting 
the computed turn, based on maximum allowable values of lateral acceleration 
and jerk. At present, these maximums are internally set to infinity and thus, 
ignored in the computation. 


Available literature on passenger comfort, however, did not provide a means of 
evaluating the interactive effect of acceleration and jerk occurring 
simultaneously. Therefore, the following relationship was assumed because of 
its simplicity: 


G + J ^ 1 

g max Jmax 

Data on passenger comfort limitations are difficult to find. Data for trains 
indicate that lateral G's of 0.12 and lateral J’s of 0.55 are comfortable to 
90 percent of the passengers and rated acceptable by 95 percent. Higher 
levels are probably acceptable on a taxiing aircraft because of the short 
duration and improved seat design. 

Equations are included to limit the rate of change in turn radius to 
acceptable levels for passenger comfort. These limits are presently 

internally set to infinity, but the equations follow the R dot equation and 
are : 


2 

n y = V /(g c R) = lateral acceleration 

J available = J max x ^ ^ G / G max^ 

R dot = ■■ g v y j / 

max g c x k x J avallable / V 

Pro bab ilistic Input 

Operation of the Model requires two data sets, named HISPEX and SND. HISPEX 
is an input data set of landing characteristics and their related deviations 
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supplied in the following format. As seen below there are only 3 cards in the 
input file, making it necessary to edit the file with every parameter change. 
The first line of 80 or less characters in a 20A4 format, is a title line to 
identify the data set. Lines two and three are landing characteristic data 
inputs, in a F5.0 format, identified as follows: 





(Any user-desired message - 

FORMAT 


TITLE CARD 


up to 80 characters) 

20A4 

UA 

US 

VA 

VS WA WS 

6F5.0 

XA 

XS 

YS 

CDIST ZA NTIP 

6F5.0 


Card 1: 

1. Title Card - any user desired message may be input, up to 80 characters in 
length. The message is printed with each run to allow the user to 
identify various cases. 

Card 2: 

1. UA - Touchdown location: Distance from runway threshold to point of 

touchdown (Feet). 

2. US - Standard deviation of touchdown location, used in calculating average 
exit speed standard deviation (Feet). 

3. VA - Touchdown speed - Average speed of aircraft between threshold and 
touchdown. Since flight test data indicates that a decrease in velocity 
from threshold to touchdown is usually on the order of 5 knots, the speed 
is averaged between the threshold and the start of deceleration (ft/sec). 

4. VS - Standard deviation of touchdown speed (ft/sec). 

5. WA - Time from touchdown to start of deceleration. Used to calculate 
speeds and occupancy time (seconds). 

6. WS - Standard deviation of WA. Accounts for discrepancy of pilot 
technique in time to apply brakes and reverse thrust (seconds). 
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Card 3: 


1. XA Deceleration rate - Rate of change of aircraft speed from start of 
deceleration to start of exit (ft/sec^). 

2. XS - Standard deviation of deceleration rate. Accounts for discrepancy of 
pilot technique in aircraft deceleration (ft/sec^). 

3. YS - Standard deviation of 50% speed used in probability calculation. A 
value of 1.0 is acceptable for most cases (ft/sec). 

4. CDIST - Exit location. Distance from threshold to point of which aircraft 
begins exit maneuver. This is the parameter to be varied to obtain a 
range of occupancy times. The solution for a specific set of input data 
is the exit location that satisfies a 40 second occupancy time and 100% 
probability (feet). 

5. ZA - Speed at which reverse thrust is shut off, which calculates B 
distance on runway. (See Figure 2, page 7) This is used to check that 
reverse thrust is terminated prior to beginning exit maneuver, as it would 
not be practical to consider a case in which the pilot must be burdened 
with a workload of reverse thrust and exit maneuvering simultaneously 
(ft/sec). 

6. NTIP - Distance from aircraft nose to wingtip, as measured along 

longitudinal axis of aircraft. (See Item 3, page 24) This parameter was 
added to the program to calculate the fraction of time the aircraft 
occupies the runway prior to wingtip over threshold. Since the optimum 
path program treats the aircraft as a single point located on its wingtip, 
occupancy time prior to wingtip occupancy is calculated by: (feet). 

Time = distance/speed => NTIP/Touchdown Speed 
Table 8 shows the 5 terps category HISPEX. DAT inputs. 
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TABLE 8 


High Speed Exit Data 
(HISPEX.DAT) 


CESSNA 172 HAWK 
500, 10, 97, 5, 3, 1 

3.5, 1, 1, 1.1777, 60, 27 


TERF’S CATEGORY B - CESSNA CITATION SERIES 500 
1000, 20, 184, 10, 4, 1 

6, 1, 1, 1.3762, 140, 19 


TERRS CATEGORY C - BOEING 747-100 
1500, 30, 235, 10, 5, 1 

6, 1, 1, 1.4075, 150, 144 


TERRS CATEGORY C - DC- 10 SERIES 10 
1500, 30, 230, 10, 5, 1 

6, 1, 1, 1.5216, 150, 130 


TERRS CATEGORY D CONCORD 
1500, 30, 275, 15, 5, 1 

6, 1, 1, 1.6300, 150, 150 


TERRS CATEGORY E - LOCKHEED F-104 
1500, 30, 370, 25, 5, 1 

6 , 2 , 1 , 1 . 8400 , 1 50 , 35 


Input data set SND is a standard table of normal distribution values. The 
program reads in Table 9, as 10 lines in a 10F5.4 format, and calculates the 
probability value, given a Z value as input. These Z values are calculated at 
various points in the model by the standard probabilistic techniques: 

Z(x) = x - x 

O 
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TABLE 9 


Standard Normal Probability Table 
(SND.DAT) 


.0000. 

. 0398 . 

. 0793 . 

. 1179 . 
. i 554. 
. 1915. 
,2257. 
.2580 . 
.2881 . 
.3159. 
, 3413 . 
. 3643 . 
.3849. 
.4032. 
.4192. 
.4332. 
.4452. 
.4554. 
.4641 . 
.4713. 
. 4772 . 
.4821 . 


0040 . 
0438. 
0832. 
1217. 
1591 . 
1950 . 
2291 . 
2611 . 
2910 . 
3186. 
3430. 
3665. 
3869. 
4049. 
4207. 
434'.:. 
4463 . 
4564. 
4649. 
4719. 
4778. 
4826. 


0080 
0478 
0871 
1 255 
1628 
1985 
2324 
2642 
2939 
3212 
3461 
3686 
3883 
4 066 
4222 
4357 
4474 
4573 
4656 
4726 


. 0 120 . 
. G51 7 . 
.0910 . 
. 1293. 
. 1664 . 
.2019 . 
. 2357 . 
.2673. 
. 2967 . 
. 3238 . 
.3485 . 
. 3708 . 
.3907 . 
.4082. 
.4236. 
.4370 . 
.4484. 
.4582. 
.4664 . 
.4732. 


0160.0199 
0557.0596. 
0948 . 0987 
1331 . 1368 
170 0 
2054 
2389 
2704.2734 
2995.3023 
3264.3289 
3508.3531 
3729 
3925 
4 0 99 


4251 
4382 
4495 
4591 

4 671 
4738 
.4783.4793 
4838 


4830 .4834 


. 0239. 

0636. 
. 1026. 
. 1406 . 
. 1736.1772. 
.2088.2123. 
.2422.2454. 
.2734.2764. 
.3051 . 
.3315. 
.3554. 
.3749.3770. 
.3944.3962. 
.4115.4131 . 
.4265.4279. 
.4394.4406. 
.4505.4515. 
.4599.4603. 
. 4673 . 4686. 
.4744.4750 . 
.4798.4803. 
.4342.4846. 


0279 
0 675 
10 64 
1443 
1808 
2157 
2486 
2794 
30 78 
3340 
3577 
3790 
3980 
4147 
4292 
4418 

Acrpc: 


. 0319 
. 0714 
.1103 
. 1480 
. 1844 
.2190 
.2517 
.2823 
.3106 
.3365 
. 3599 
.3810 
.3997 
.4162 
. 4306 
. 4429 
. 4535 
4616.4625 
4693.4699 
4756.4761 
4808 
4350 


. 4812 
.4354 


. 0359 
. 0753 
. 1141 
. 1517 
. 1879 
. 2224 
.2549 
. 2852 
.3133 
. 3389 
.3621 
.3830 
.4015 
.4177 
.4319 
. 4441 
. 4545 
.4633 
.4706 
.4767 
4317 
4857 


. 4861 . 4864 . 4860 . 4871 . 4875 . 4878 . 4831 . 4834 . 48B7 . 4890 
.4893.4896.4898.4901 .4904.4906.4909.4911 .4913.4916 
.4918 . 4920 . 4922 . 4925 . 4927 . 4929 . 4931 . 4932 . 4934 . 4936 
. 4938 . 4940 . 4941 . 4943 . 4945 . 4946 . 4943 . 4949 . 4951 . 4952 
. 4953 . 4955 . 4956 . 4957 . 4959 . 4960 . 4961 . 4962 . 4963 . 4964 
.4965 .4966 . 4967 . 4963 . 4969 . 4970 . 4971 .4972 . 4973 . 4974 
. 4974 . 4975 . 4976 . 4977 . 4977 . 4978 . 4979 . 4979 . 4980 . 4931 


. 4981 . 4982 . 4902 . 4983 . 4984 . 4984 . 4985 . 4905 . 4986.4986 
. 4987 . 4987 . 4937 . 4988 . 4983 . 4989 . 4989 . 4989 . 499 0 .4990 
. 4990 . 4991 . 4991 . 4991 . 4992 . 4992 . 4992 . 4992 . 4993 . 4993 
. 4993 . 4993 . 4994 . 4994 . 4994 . 4994 . 4994 . 4995 . 4995 . 4995 
. 4995 . 4995 . 4995 . 4996 . 4996 . 4996 . 4996 . 4996 . 4996 . 4997 
. 4997 . 4997 . 4997 . 4997 . 4997 . 4997 . 4997 . 4997 . 4997 . 4998 
. 4998 . 4998 . 4993 . 4998 . 4998 . 4998 . 4998 . 4998 . 4998 . 4998 
. 4998 . 4998 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 
. 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 
. 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 . 4999 
. 5000 . 5000 . 5000 . 5000 . 5000 . 5000 . 5000 . 5000 . 5000.5000 
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Incrementing by one, through the discrete probability and summing, minimum 
occupancy time, and arbitrary occupancy time, the program iteratively 
calculates values until the final value conditions are met. 


Exit Path Subroutine Input 

The subroutine named EXPATH tracks the aircraft through the exiting maneuver. 
The only input required by EXPATH that is dependent on the probabilistic model 
is exit speed (speed @ point C). A separate data file, named (PATH.DAT) is 
input to the path model, consisting of the data in Table 10. 


TABLE 10 

INPUT DATA TO EXIT PATH COMPUTATION 


A/C Weight 

Yaw Inertia 

Wheelbase 

% Weight on 
Main Gear 

Theta Max 

YMU 

Deceleration 


Runway Width 

Wingspan 

Distance to 
Taxi way 


Y MU Scrub (1) 

Y MU Scrub (N) 

Radius (1) 
Radius (N) 


F5. 


FORMAT 

F12.3 


F12.3 


F12.3 


, G10.1 


The input data required is described below: 


1. Aircraft weight - Should be maximum landing weight to obtain critical case 
(pounds). 

2. YAW Inertia - Resistance to movement about vertical axis of aircraft 
(slug-Ft 2 ). 
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3. Wheelbase - Distance from nose gear to main gear, along longitudinal axis 
of aircraft (feet). 

4. Percent weight on main gear - Percent of total aircraft weight on main 
gear, in maximum landing weight configuration. 

5. THETA MAX - The angle at which the wingtip and stabilizer tip are both 
touching the runway edge. Figure 4. The angles for all 5 Terps category 
aircraft are shown in Table 11. The tracking program checks to see that 
this angle is not exceeded during the exit maneuver. If the angle is 
exceeded, the program is terminated and an appropriate output message is 
printed on the screen. 



FIGURE 4. SKETCH OF THETA MAX 


6. Y MU MAX - Measure of maximum resistance offered by nose gear. A MU (Cng) 
of 0.2 is considered to be a practical limit for maneuvering at speeds of 
60 - 90 feet per second. 

7. Deceleration rate - This is the deceleration rate employed in the model 
during the turning maneuver. This variable allows the analyst to vary 
deceleration on the runway and through the turn (ft/sec 2 ). 
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Tabic 11 


Exit Path Data 
(PATH. DAT) 


CESSNA 172 HAWK 


2550. 2308 . 

41. 0.2 

150. 35.9 

0.00 i.CE^i 

0.02 50 . 


5. A 
- 2.0 
600 . 


75 . 


TERPS CATEGORY B - CESSNA CITATION SERIES 500 


11000. 59489 . 15.7 88. 

38. 0.2 -2.0 

150. 43.8 600. 

0.00 i.OE+31 

0.02 50. 


TERPS CATEGORY C - BOEING 747-100 

564000. 42000000. 84. 96.4 

45. 0.2 -2.0 

150. 195.7 600. 


00 

0 . 1EH32 

01 

5287. 

019 

2638. 

039 

1318 , 

059 

877. 

079 

657. 

099 

524. 

153 

347. 

212 

257. 

353 

167. 

548 

121 . 

793 

99. 


TERPS CATEGORY C - DC-10 SERIES 10 


363500. 16350338. 
44. 0.2 

150. 155.3 


0 .00 

1 .0EH31 

0.01 

5000 . 

0 . 02 

2000 . 

0 . 025 

1000 . 

0 . 027 

800 . 

0.03 

200 . 


72.5 93.4 

- 2.0 

600 . 
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Tabic 11 (Cont. ) 

TERPS CATEGORY D - CONCORD 


2450 0 0 . 12664S09. 

37 . 0.2 

ISO. 03.0 

0.0 0 0 . i E-i 32 

0.01 i 0 0 0 . 

0 . 02 80 0 . 

0 . 03 50 . 


TERRS CATEGORY E - LOCKHEED F-104 


22000. 137776. 

20 . 0.2 

150. 21.9 

0.00 O.iE+32 

0.02 50. 


8. Runway width - The width of the runway surface (feet). 

9. Wingspan - Total straight line distance from aircraft wingtip to wingtip 
(feet) . 

10. Distance to Taxiway - The distance between the centerlines of the runway 
and a parallel taxiway (Feet). 


11. YUU Scrub and Radius - These values are the scrub function as determined 
by Fig. 3. They vary for different aircraft type but are negligible for 
small aircraft with a dual main gear configuration. 


Clearing Distance 


Variables 8 and 9 are entered into the program to determine when the 
aircraft is clear of the runway. During the tracking calculation, the 
aircraft is regarded as a single point located on its wingtip. From 
Figure 5, it can be seen that as an aircraft proceeds down the runway, the 


15.0 92.3 

- 2.0 
600 . 


59.7 94 . 7 

- 2.0 
600 . 
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RUNWAY EDGE 



total distance necessary for the wingtip to travel in order to clear the 
runway can be given by: 

clearing distance = ( Runway width + Wingspan ) 

2 


At intervals of 1/1000 of a second, the model calculates the distance the 
wingtip has traveled laterally and compares this value to (RW + WS)/2. 

When Y > (RW + WS)/2, the model checks to see if Theta < Theta max, then 
returns the time lapsed to the probabilistic model (feet). 
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Output symbols for the Probabilistic Model 


1. ASCSV - speed at the start of the exit maneuver (ft/sec). 

2. ACS - speed at the end of the exit maneuver (ft/sec). 

3. PTIME - time on exit (curved) path, to clearing the runway surface 
(seconds). 

4. CDIST - distance from threshold to the beginning of the exit maneuver 
(feet) . 

5. XDIST - distance traveled during the exit maneuver, to clearing the runway 
surface only (feet). 

6. TDIST - total distance traveled, CDIST plus XDIST (feet). 

7. PROB - probability of using the high speed exit (percent). 

8. AVG TIME - average time summed through the probability step incrementation 
(seconds). 

9. STD DEV — standard deviation of calculated occupancy times. 

10. XTIME - time from the threshold to the beginning of exit path (seconds). 

11. TIME - total occupancy time, from threshold to clearing the runway surface 
(seconds) . 

12. YA - probability speed (ft/sec). 

Output symbols for exit path 

1. Time - time in seconds from the start of the exit path. 

X ~ X— coordinate distance along the longitudinal axis of the runway from 
the beginning of the exit path (feet), Fig 6. 
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( 0 . 0 ) 


X 


<£ RUNWAY 



FIGURE 6. EXIT PATH COORDINATES 


3* Y - Y-coordinate distance along the lateral axis of the runway from the 
beginning of the exit path (feet). 

4. THETA - angle produced between X and Y coordinate (degree). 

5. Radius - instantaneous radius of curve which would fit between two 
consecutive points along the path (feet). 

6. R rate - instantaneous change in radius calculations. 

7. Mu Scrub - coefficient of friction contribution from tire scrub - C _ 

s * 

nose gear. 

8. Mu Cent - coefficient of friction contribution from centrifugal forces 
(percent) - C a y, nose gear. 

9. Mu I - coefficient of friction contribution from rotational inertia of the 
aircraft - C^.. 

10. G cent - centrifugal G forces. 

11. Velocity - instantaneous speed at specific time (ft. sec). 
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RUNWAY TURNOFF MODEL APPLICATIONS 


The following examples consider five Terps categories, as example aircraft 
which were used to run the program based on each of these aircraft's 
parameters. 

Terps Category A 

1. Cessna 172 Hawk 

Terps Category B 

1. Cessna Citation series 500 

Terps Category C 

1. Douglas DC-10 series 10 

2. Boeing 747 series 100 

Terps Category D 

1. British/French Concord 

Terps Category E 

1. Lockheed F-104 

Tables 12-23 are produced from the input data specified in Tables 8, 9, and 
11. Speeds, times, probabilities, distances, and exit paths are shown for 
each following example. 

A complete listing of the Fortran Programming can be found in the Appendix. 
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Table 12 

TERPS CATEGORY A 
CESSNA 172 HAWK 

NASA HIGH SPEED EXIT ANALYSIS 
INPUT SUMMARY 

PROBABILISTIC MODEL 


TOUCHDOWN LOCATION 
TOUCHDOWN SPEED 
TIME-TOUCHDOWN TO DECEL 
DECELERATION RATE 
507. PROBABILITY SPEED 
DISTs THRESHOLD-START OF TURN 
SPEED AT MAJOR DECEL REDUCE 
NOSE-WINGTIP DIST 


AVERAGE 

STD DEV 

UNITS 

500. 

10. 

FEET 

97. 

5. 

FT/SEC 

3. 

1 . 

SECONDS 

3.50 

• 

o 

o 

FT/SEC**2 

65. 

1 . 

FT/SEC 

1777. 


FEET 

60. 


FT/SEC 

27. 


FEET 


AIRCRAFT TRACKING PROGRAM 


AIRCRAFT CHARACTERISTICS 
WEIGHT (LBS) YAW 

2550.0 


INERTIA 

2308.0 


WHEELBASE 

5.4 


X LOAD ON MAIN GEAR 
75.0 


FLIGHT CONDITIONAL 
EXIT VELOCITY 
50.070 


PARAMETERS 
MAX EXIT ANGLE 
41.000 


MAX Y MU DECEL RATE 

0.200 - 2.000 


SYSTEM WIDTHS 

RUNWAY WINGSPAN 

150.00 35.90 


CLEARING DIST. 
92.95 


DIST TO TAXIWAY 
600.00 


SCRUB CHARACTERISTICS 
Y MU RADIUS 

0.0000 O. 1E+32 
0.0200 50.0 
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TABLE 13 


iLCUl ATl!l» RESULTS 


CRITICAL VALUES i 
ASC5V 

ASC 

PTIHE 

CDXST 

50.07 

36.51 

6.79 

1777.00 


XDIST TDIST 

26V. 04 2046.04 


PROB 

AVC TIME 

STD DEV 

XTInE 

TTIME 

YA 

1.0000 

21.5632 

2 . 9 0 C 4 

32.8747 

39.6604 

65.000044 

O.V906 

21.8324 

2.9004 

33.143? 

39.9296 

64 .000044 


CONSTANT NLC SIDE LOAD ANALYSIS 
PASSENGER COMFORT ICNOnCP 


TIME 

X 

Y 

THETA 

RADIUS 

R RATE 

MU 

SCRUB 

MU 

CENT 

0 

.000 

0 , 

.000 

0.000 

0 .1300 

100000.00“ 

40041660 

.00 

0 

. 020 

0 

.001 

0 

t; 

.500 

;ooo 

24 

* ' 4 ?: 

.704 
.048 ’ 

0.173 
” " 1. 126 " 

1 .122 
”3.521 * 

724 .94 
500 . 4 if 

-920 

“213' 

.69 
.'4 0 

0 

“‘O' 

. 020 
. 020 ' ’ 

0 

O' 

.103 
.14 3 

i 

.500 

72 

.741 

3. 183 

6.482 

431 .20 

-90 

.23 

0 

.020 

0 

.160 

2 

.000 

95 

.793 

6.453 

9.717 

396.92 

-53 

.26 

0 

.020 

0 

. 166 

2, 

,S00 

110 

. 126 

10.954 

13.108 

374,21 

-3? 

.63 

0 

.020 

0 

.169 

3 

,000 

139 

.665 

16.661 

16 . tno 

356,01 

-33 

.95 

0 

.020 

0 

.169 

3, 

,500 

160 , 

.336 

23.531 

20 .198 

339.70 

-31 

.25 

0 

.020 

o 

. 170 

4. 

.000 

1G0 , 

, C 67 

31 .50? 

23.870 

324.57 

-2? 

.73 

0 

.020 

0 

.169 

4. 

.500 

190, 

,790 

40.532 

27.625 

309.97 

-28 

.69 

0 

.020 

0, 

.169 

5, 

.000 

216 , 

.439 

50.527 

31.463 

29S.04 

-27 

.86 

0 

. 020 

0 , 

,169 

5. 

.500 

232, 

,954 

61.41? 

35.306 

282.10 

-27 

.11 

0 

.020 

0, 

.168 

6. 

.000 

240, 

.279 

73.125 

3V .4 03 

268.72 

-26 

.40 

0 , 

.020 

0, 

.168 

6, 

.500 

262 . 

.362 

85.554 

43.506 

255.69 

-25 

.70 

0, 

. 020 

0 , 

,167 

7, 

.000 

275, 

. 158 

98.614 

47.7C9 

243.01 

-25 

.02 

0, 

.020 

0 , 

,166 

7, 

,500 

286, 

.620 

112.206 

52.013 

230.67 

-24 

.33 

0 

.020 

0, 

,166 

0. 

,000 

296 , 

.741 

126.223 

56.422 

218.68 

-23 

.65 

0, 

. 020 

0. 

. 165 

0, 

,500 

3G5 

.471 

140.55? 

61.942 

207.02 

-22, 

.97 

0 

.020 

0 , 

, 164 

7, 

,030 

312, 

.805 

155.099 

65 .678 

195.70 

-22 

.29 

0 , 

. 020 

0. 

,163 

9. 

.500 

3 1 B . 

.736 

169.726 

70. 335 

184.72 

-21 , 

,63 

0 , 

. 020 

0. 

, 162 

10 , 

.000 

323, 

,269 

104.323 

75.21? 

174.00 

-20 , 

.95 

0. 

, 020 

0 . 

,161 

10, 

,501 

326, 

,420 

198.767 

00.237 

163.77 

-20, 

.20 

0 , 

.020 

0 . 

, 160 

11. 

.001 

323, 

,214 

212,938 

05.394 

1S3.80 

-19, 

,61 

0. 

. 020 

0 . 

159 

11 . 

.501 

320, 

,692 

226.714 . 

90.690 

144.16 

-10, 

,93 

0 . 

,020 

0 . 

158 

12, 

,001 

327. 

.906 

239.974 

96.157 

134.86 

-18. 

,26 

0 . 

020 

0. 

157 

12, 

.501 

325, 

,922 

252.602 

101.778 

125.90 

-17. 

,59 

0. 

020 

0. 

155 

13, 

,002 

322. 

,819 

264,406 

107.571 

.117.27 

-16. 

92 

0 . 

020 

0. 

154 

13. 

,502 

310. 

,690 

275.521 

113,543 

1 08, ?B 

-16. 

2*6 

0 . 

020 

0. 

152 

14. 

,002 

313. 

642 

205.611 

119.706 

101.01 

-IS. 

60 

0 . 

020 

0 . 

ISO 

14. 

, S02 

307 . 

795 

294.670 

126. C67 

93.33 

-14. 

94 

0 . 

020 

0 . 

148 

15. 

,002 

301 . 

202 

302.626 

132.639 

86.07 

-14. 

29 

0 . 

020 

0. 

145 

15. 

,503 

294 . 

,24? 

309.423 

139.433 

79.10 

-13. 

63 

0. 

020 

0. 

143 

16. 

,003 

286. 

,848 

315.023 

146.459 

72.44 

-12. 

90 

0. 

020 

0 . 

140 

16. 

,503 

279. 

,241 

319.40? 

153.729 

66.12 

-12. 

33 

0 . 

020 

0 . 

137 

17. 

,002 

271 . 

594 

322.587 

161.257 

60.11 

-11 . 

70 

0. 

020 

0 . 

133 

17. 

502 

264 . 

076 

324.590 

169.052 

54.42 

-11 . 

06 

0 . 

020 

0 . 

130 

4440UT OF 

TABLE 

BOUNDS444 









^ I C CENT VELOCITY 


0.179 0. 

0.077 0. 

“0.037 * 0. 

0.020 0 . 

0.014 0. 

0.011 0 . 

0.011 0 . 

0.010 o. 

0.011 0 . 

0.011 o. 

0.011 0 . 

0.012 0 . 

0.012 0 . 

0.013 0. 

0.014 0. 

0.014 0. 

0.01S 0. 

0.016 0. 

0.017 0. 

0.018 0. 

0.019 0. 

0.020 C. 

0.021 0 . 

0.022 0 . 

0.023 0. 

0.C2S 0, 

0.026 0. 

0.02G 0. 

0.030 0. 

0.032 0. 

0.03S 0. 

0.037 0. 

0.040 0. 

0.043 0. 

0.047 0. 

0.050 0. 


001 50.070 
103 49.070 
143 48.071 
160 47.072 

166 46.072 
169 45.073 

169 44.073 

170 43.074 
169 42.074 
169 41.07S 
169 40.075 
160 39.076 

160 33.077 

167 37. 77 
166 36.078 
166 35.078 
165 34.077 
164 33.07? 
163 32. *80 
162 31.080 

161 30.07? 
160 29.079 
159 2B.076 
1SB 27. 076 
157 26.077 
155 2S.077 
154 24.077 
152 23.076 
150 22.076 
140 21.075 
145 20.075 
143 19.075 
140 18.074 
137 17.074 
133 16.073 
130 15.073 
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TABLE 14 


TERPS CATEGORY B - CESSNA CITATION SERIES 500 

NASA HIGH SPEED EXIT ANALYSIS 
INPUT SUMMARY 

PROBABILISTIC MODEL 

AVERAGE STD DEV UNITS 

TOUCHDOWN LOCATION 
TOUCHDOWN SPEED 
TIME-TOUCHDOWN TO DECEL 
DECELERATION RATE 
507. PROBABILITY SPEED 
DIST: THRESHOLD-START OF TURN 
SPEED AT MAJOR DECEL REDUCE 
NOSE-WINGTIP DIST 


AIRCRAFT TRACKING PROGRAM 


AIRCRAFT CHARACTERISTICS 

WEIGHT (LBS) YAW INERTIA 

11000.0 59489.0 

WHEELBASE 

15.7 

7. LOAD ON MAIN GEAR 
88.0 

FLIGHT CONDITIONAL 
EXIT VELOCITY 
97.693 

PARAMETSRS 

MAX EXIT ANGLE 
38-000 

MAX Y MU 
0.200 

DECEL RATE 
-2.000 

SYSTEM WIDTHS 
RUNWAY 
150.00 

WINGSPAN 

43.80 

CLEARING DIST. 
96.90 

DIST TO TAXIWAY 
600.00 


SCRUB CHARACTERISTICS 
Y MU RADIUS 


0.0000 0.1E+32 

0.0200 50.0 


1000. 

20. 

FEET 

184. 

10. 

FT/SEC 

4. 

1 . 

SECONDS 

6.00 

1.00 

FT/SEC0ND**2 

113. 

1 . 

FT/SEC 

3762. 


FEET 

140. 


FT/SEC 

19. 


FEET 
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TABLE 15 


cfA cuL/urn Reruns 


CRITICAL VALUfS: 


AS CSV 

V7.6V 

asc 

03.65 


Mire 

7. 0 S 

rm:i 
JVM*. 0 0 

VDIT.T 
t. 4 

1 S 1 V f 
< ;*.t . 

PROS 

A VC T1HE S 

rTE> rev 

XT IPX 

TTIME 

YA 


i .0000 
0.999? 

23.B997 

23.0??7 

2.324? 

2.324V 

32 . 966? 
32.966? 

37.9935 

39.9935 

113.00004* 
1 12.0000** 



CONSTANT NLC SIDE 

LOAD ANALYS 

IS 








PASSENCEP COMKORT 

IGNORED 









TIME 

X 

Y 

THETA 

RADIUS 

R RATE 

MU SCRUB 

MU CENT 

MU I 

C CENT 

t'ELO: 

0 .000 

0 .000 

0.000 

0.000 

103000.00 

“555541 0 . 00 

0.020 

0.003 

0.177 

0.003 

97.693 

0 .500 

48.596 

0.110 

0.364 

4256.54 

-6421 .15 

0.020 

0 . 063 

0.112 

0 .068 

96.694 

1.000 

96.688 

0.740 

1.224 

2637.92 

-1609.30 

0.020 

0.103 

0 . 072 

0 . 1CB 

95.694 

1.500 

144.261 

2 . 232 

2 . 375 

2100.55 

-689.03 

0 . 020 

0 . 132 

0 .048 

0.132 

94.695 

2.000 

191.290 

4.751 

3 . 7o8 

1054.76 

-371 .53 

0.020 

0 . 147 

0.033 

0.147 

93.695 

2.50 0 

237.742 

8.413 

5 272 

1700.57 

-229.93 

0.020 

0.1S6 

0.024 

0.156 

92.696 

3.0C0 

283.580 

13.27? 

6.B65 

1613.43 

-157.99 

0.020 

0 . 162 

0.010 

0 . 162 

91 .697 

3.500 

328.767 

19.377 

e .521 

1545.23 

-118.44 

0 . 020 

0.165 

0.015 

0 . US 

90 .697 

4.000 

373.263 

26.716 

10.224 

1492.22 

-95.54 

0.020 

0.167 

0.013 

0.167 

D9.698 

4.500 

417.027 

35.292 

11.^62 

14*38. 16 

-Bi .79 

0.020 

0 .169 

0.011 

0.169 

89.690 

5.000 

460. C21 

45.092 

13.731 

1409.56 

-73.26 

0.020 

0.169 

0 .011 

0.169 

87.699 

5.500 

502.206 

56.099 

15.526 

1374.39 

-67.79 

0 . 020 

0.170 

0.010 

0 . 170 

85.699 

6.003 

543.544 

60.288 

17.345 

1341.46 

-64 .17 

0.020 

0.170 

0 .010 

0 .170 

85.700 

6.50C 

553.997 

61 .635 

1° . 1P6 

1310.03 

-61 .67 

0.020 

0.170 

C.010 

0 .170 

84.701 

•*.009 

623.526 

76.111 

21 . 04? 

1279.60 

-59.04 

o.oro 

0.170 

0.013 

C.170 

83 v ; 1 

.500 

662.101 

111 .684 

22.733 

1250 .12 

-5e.45 

0.020 

0.170 

0.010 

C.170 

62 .702 

8. 000 

699.682 

120.323 

24 .839 

1221 . 19 

-57.31 

0.020 

0 .170 

0.010 

0 .170 

81 .702 

e.soo 

736.234 

145.991 

2 1> , V 6 7 

1192.78 

-56.34 

3,020 

0.170 

0.010 

0 . 17C 

60 .703 

9.000 

771.727 

164.652 

23.716 

1164.83 

-55.40 

0.020 

C .169 

0.011 

0.16? 

79.703 

9.500 

806.125 

184.270 

33.667 

1137.29 

-54.69 

0.020 

0.169 

o 

o 

0 .169 

7B.704 

10.000 

839.400 

204 . 803 

32 . 6c 1 

1110.14 

-53.92 

0.020 

0.169 

0.011 

0.169 

77.704 

10.501 

871.51? 

226.210 

34.697 

1083.37 

-53.19 

0 . 020 

0.1 6 9 

0.011 

0.169 

76. 70S 

11.001 

902.455 

240.449 

36.737 

1056.95 

-52.47 

0 . 020 

0.163 

0.012 

0 .168 

75.7:5 

11.501 

932.100 

271.475 

39,83) 

1030.09 

-51.77 

0 . D20 

0 .168 

0 .012 

0 . 168 

74.706 

12.001 

960.668 

295.243 

4 C . BVO 

1005.10 

-51 . 07 

0 . 020 

0 .163 

0.012 

0.168 

73.707 

12.501 

987.094 

319.705 

43.CC3 

979.82 

-50.30 

0.020 

0.168 

0 .01? 

0.168 

72.707 

13.002 

1013.834 

344.013 

45.141 

954.03 

-49.69 

0.020 

0.167 

0.013 

0 . 167 

71 .708 

13.502 

1038.468 

370.516 

4 7 .3 36 

930 .12 

-49.01 

0.020 

0 .167 

0.013 

0.167 

70.706 

14.002 

1061 .775 

396.762 

49.497 

905.79 

-40.32 

0 . 020 

C .167 

0.013 

0.167 

69.70? 

14.502 

1003.739 

423. 50U 

51.715 

e91 .80 

-47.64 

0.020 

0 . 166 

0.014 

0.166 

68.70? 

IS. 002 

1104.342 

450 .674 

53.961 

058.14 

-46 . 96 

0.020 

0 . 166 

0 .014 

0.166 

67.710 

15.503 

1123.572 

473.230 

56.236 

834.04 

-46.23 

0.020 

0 .166 

0.014 

0 .166 

66.710 

16.003 

1141.416 

506.112 

50.540 

811 .07 

-45.60 

0 . 020 

0 . 16S 

0.015 

0 . 165 

65.711 

16.503 

1157.865 

534.261 

60 .e/3 

709.24 

-44.92 

0.020 

0.165 

0.015 

0.165 

64.712 

17.002 

1172.914 

562.619 

63.228 

766.95 

-44 .24 

0 . 020 

0 .164 

0 .016 

0 . 164 

63.712 

17.502 

1186.555 

591 .126 

65.633 

745.00 

-43.56 

0 . 020 

0 .164 

0.016 

0.164 

62.713 

FINAL 

VALUES 










17.658 

1190.524 

600.042 

65.633 

738.22 

-43.35 

0.020 

0.164 

0.016 

0.164 

62.401 


:ty 
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TABLE 16 


TERPS CATEGORY C - DC-10 SERIES 10 

NASA HIGH SPEED EXIT ANALYSIS 
INPUT SUMMARY 

PROBABILISTIC MODEL. 


TOUCHDOWN LOCATION 
TOUCHDOWN SPEED 
TIME-TOUCHDOWN TO DECEL. 
DECELERATION RATE 
50% PROBABILITY SPEED 
DIST: THRESHOLD-START OF TURN 
SPEED AT MAJOR DECEL REDUCE 
NOSE-UINGTIP DIST 


AVERAGE 

STD DEV 

UNITS 

1500 . 

30 . 

FEET 

230 , 

10 . 

FT/SEC 

5 * 

1 . 

SECONDS 

6.00 

1 .00 

FT/SEC*: 

164. 

1 . 

FT/SEC 

5216. 


FEET 

150 . 


FT/SEC 

130 . 


FEET 


AIRCRAFT TRACKING PROGRAM 


AIRCRAFT CHARACTERISTICS 


WEIGHT ( LBS ) 
363500 . 0 


YAW INERTIA 
16350338.0 


WHEELBASE 

72.5 


% LOAD ON MAIN GE^tR 
93./. 


FLIGHT CONDITIONAL 
EXIT VELOCITY 
148.688 


PARAMETERS 
MAX EXIT ANGLE 
44.000 


MAX Y MU 

0.200 


DECEL RATE 

- 2.000 


SYSTEM WIDTHS 

RUNWAY 
150 .00 


WINGSPAN CLEAR DIST. 

155.30 152.65 


DIST TO TAXI WAY 
600.00 


SCRUB CHARACTERISTICS 


Y MU 

RADIUS 

0.0000 

0 . 1E+32 

0.0100 

5000 . 0 

0 .0200 

2000.0 

0.0250 

1000.0 

0 . 0270 

800 . 0 

0 . 0300 

200.0 
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TABLE 17 


tAt UJl All n Kl SIM l!i 
J T 1 CAL VAUILSl 


ASL'JV 

A!;r. 

p r i r.r; 

mm 

XI 1ST 

ibi'.r 

140.6? 

130.30 

9 . 20 

\»:*t a.uo 

1 26U . 66 

64U4 .66 


PROP 

AVC TIHE 

sto Drv 

XTIHE 

TT 1HE 

YA 

1.0000 

0.99?? 

2b. 6JV7 

i . 3 pv; 

1 . 3255 

to . :»•> a 

3Ct .ltUV4 

40. on?5 
40.0O7S 

Jf i. v ** o »• i * 
1 /. i . 0 0 t 0 1 * 


CONSTANT NLC SIDE 
PASSENCER COr.roRT 

TIHE X 

LOAD ANALYSIS 
IGNORED 

Y THE f A 

RADIUS 

R RATE 

HU SCRUB 

HU CENT 

HU I 

C CENT 

UELCT.ITY 

0.000 

0 .000 

0.000 

0 .003 

l oopco .oo 

-1223773.37 

0.010 

0.007 

0 .103 

0.007 

1 48 ,6Df 

o .soo 

74 . 093 

0 .000 

C .163 

15359 ,7V 

>23157. 03 

0.010 

0 . 044 

0 .146 

0 . 044 

1 47.688 

1.000 

147.6HS 

. 0.SC8 

0 .534 

9143.11 

-6629.52 

0.010 

0.073 

0.117 

0.073 

14C.60? 

1 .500 

220.773 

i .514 

1 . 069 

6092.57 

-3057.40 

0.010 

0.076 

0 . 074 

0 . 096 

145.609 

2.000 

293.345 

3.277 

1.735 

5739 . 13 

-173S. 16 

0 .010 

0.113 

0 .077 

0.113 

144.690 

2.500 

365.390 

5.934 

2.5C6 

5045.25 

-1107.27 

0.010 

0.127 

0.063 

0.127 

143. 6VC 

3.000 

436.871 

9.508 

3.361 

4509.08 

-746.22 

0.011 

0.138 

0.051 

0 .138 

142.691 

3.50 0 

S07.G2Q 

14.318 

4.282 

4274.04 

-534 .10 

0.012 

0 .146 

0.042 

0.146 

141 .691 

4.000 

570.170 

20.180 

5.235 

4042.59 

-4 01 .'it 

0.013 

0.152 

0 . 035 

0 . 1S2 

1 4 ; .692 

4.50C 

647.91? 

27.214 

6.273 

3365.49 

-312.03 

0.014 

0.157 

0 . 02? 

0.157 

139.693 

5.000 

717.025 

35.440 

7.324 

3725.07 

-252.31 

0.014 

0.160 

0 . 025 

0 . 160 

138.693 

5.500 

705.471 

44.097 

0.404 

3610.20 

-209.47 

0.015 

0 .163 

0.022 

0.163 

137.674 

6.000 

853.231 

55.571 

9 .539 

3513.60 

-178. 10 

0.015 • 

0.165 

0 .020 

0.165 

126.694 

6.500 

920.278 

67.474 

1C .622 

3430.33 

-155.63 

0.015 

0.167 

0.01P 

0 .167 

1 35 . 695 

7.000 

906. S07 

80.603 

11.774 

3356.94 

-138.65 

0.01S 

0.168 

0.017 

0.168 

134.695 

7.500 

1052.131 

74.952 

12.930 

32°0 . 97 

-125.70 

0.016 

0 . 1 69 

0.016 

0 . 16? 

1 33 . 696 

B. 000 

1116.824 

110.513 

14 . 1C1 

3230.65 

-115.94 

0.016 

0.169 

0 .OiS 

0.167 

132 .696 

8.500 

11B0 .020 

127.273 

15.283 

3174.66 

-100.34 

0.016 

0 . 170 

0 .014 

0.170 

131.697 

9.000 

1243. 91S 

145.219 

16.477 

3122.03 

-102.41 

0.016 

0 .170 

0.014 

0 . 170 

120.697 

9.500 

1306.142 

164. 33S 

17..' 82 

3C72 . 03 

-97.75 

0 .016 

0 . 170 

0.014 

0 . 170 

129.698 

.000 

1367.477 

184.603 

1C .076 

3024.12 

-94.03 

0.017 

0.170 

0.013 

0.170 

120.699 

. .501 

1427.096 

206.005 

20 .120 

2977.87 

-91 .04 

0.017 

0.170 

0 . 013 

0 . 170 

127.699 

11.001 

1487,376 

228,520 

21.353 

2932 .90 

-00.61 

0.017 

0 .170 

0.013 

Q. 170 

. 126.700 

11.501 

1545.891 

252.127 

22.575 

2089.20 

-86.60 

0.017 

0 . 170 

0.013 

0.170 

125.700 

12.001 

1603.420 

276.005 

23.345 

2046.33 

-84.91 

0.017 

0 .170 

0.013 

0.170 

124.701 

12.501 

1659.940 

302.529 

25.105 

2004 . 24 

-83.40 

0.017 

0 .169 

0.013 

0 . 169 

123.701 

13.002 

1715.429 

329.277 

26.373 

2762.02 

-82 . 24 

0.017 

0.169 

0.013 

0.16? 

122.702 

13.502 

1769.065 

357. 023 

27.649 

2721 .97 

—01.16 

0 .010 

0 .169 

0 .013 

0 .169 

121.702 

14.002 

1823.227 

385.743 

20.934 

2681 .64 

-BO. 19 

0.018 

0.169 

0.014 

0 . 169 

120 .703 

14.502 

1875.475 

415.409 

30.223 

2641 .77 

- 7*? . 32 

0.018 

0.160 

0.014 

0.168 

119.704 

15.002 

1926.649 

445.995 

31 .531 

2602.31 

-72. 5i 

0.018 

0.163 

0.014 

0 .168 

1 13 704 

15.503 

1976.668 

477.473 

32.042 

2563.24 

-77 . 7 6 

0 .013 

0.163 

0.014 

0 . 163 

1 17.705 

16.003 

2025. 53S 

509.816 

34.162 

2524.55 

-77.05 

0.010 

0.168 

0.014 

0.168 

116. 70S 

16.503 

2073.231 

542.994 

35.470 

24e6. 19 

-76.37 

0.018 

0.167 

0.014 

0.167 

115.706 

17.002 

2119.739 

576.973 

36.828 

2448.17 

-75.73 

0 .019 

0.167 

O.OiS 

0.167 

114 706 

FINAL 

17.335 

VALUES 

215TK047 

600.044 

36.B28 

2423.02 

-75.30 

0.019 

0.167 

0.015 

0.167 

114.041 
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TABLE 18 


TERRS CATEGORY C - DOE I NR 747-100 

NASA HIGH SREEI) [EXIT ANALYSIS 
INPUT SUMMARY 


PROBABILISTIC MODEL 



AVERAGE 

STD DEV 

UNITS 

TOUCHDOWN LOCATION 

1500 . 

30. 

FEET 

TOUCHDOWN SPEED 

235. 

10 . 

FT/SEC 

TIME-TOUCHDOWN TO DECEL 

5. 

1 . 

SECONDS 

DECELERATION RATE 

6.00 

1.00 

ft/sec**;l 

50% PROBABILITY SPEED 
DIST: THRESHOLD-START OF TURN 
SPEED AT MAJOR DECEL REDUCE 
NOSE-WINGTIP DIST 

176. 
5125, 
150 . 
144 . 

1 . 

FT/SEC 

FEET 

FT/SEC 

FEET 


AIRCRAFT TRACKING PROGRAM 


AIRCRAFT CHARACTERISTICS 

WEIGHT (LBS) YAW INERTIA 

564090.0 42000000.0 

WHEELBASE 
84. 0 

% LOAD ON MAIN GEAR 
96.4 

FLIGHT CONDITIONAL 
EXIT VELOCITY 
160.702 

PARAMETERS 
MAX EXIT ANGLE 
43 .000 

MAX Y MU 
0.20 0 

DECEL RATE 
-2.000 

SYSTEM WIDTHS 

RUNWAY 

150.00 

WINGSPAN 

195.70 

CLEAR DIST. 
172.85 

DIST TO TAX I LAY 
600.00 


SCRUB CHARACTERISTICS 
Y MU RADIUS 


0.0000 0.1E+32 
0.0100 S287.0 
0.0190 2630.0 
0.0390 1318.0 
0.0590 877.0 
0.0790 657.0 
0.0990 524.0 
0.1530 347.0 
0.2120 257.0 
0.3530 167.0 
0.5480 121.0 
0.7930 99.0 
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TABLE 19 


CALCULATED ktSULTS 










critical vmuris: 










A5CSV 

A SC 

PTINf 

CDIST 

>:pi m 

Tt>) 

ST 



160 .70 

137.00 

n 

.46 5 

125. Oil 

16V5.1V 

65. -0. 

22 



PROB 

A VS TIME STD 

DEV 

XTIHE 

TTIKE YA 






i .0000 

24.3797 1. 

0716 

20.5S91 

40.01S1 176 

.0000** 





0.999? 

24.3797 1. 

0716 

23.5591 

40.0151 175 

. 0000** 





CONSTANT 

NLC SIDE LOAD ANALY 

SIS 








PASSENGER COMPORT 

IGNORED 


* 







TIKE 

X 

Y 

THETA 

RADIUS 

R RATE 

KU SCRUB 

KU CENT 

KU I 

C CENT 

VELOCITY 

0.000 

o.ooo 

0.000 

0 .COO 

100000 .00 

-443293.94* 

0.010 

0 .008 

0.182 

0.000 

16C.702 

0 .500 

GO. ICO 

0.055 

0.095 

32000.69 

-415C1 .80 

0.010 

0.02S 

0.163 

0.025 

159.702 

1.000 

159.701 

0.307 

0.252 

19744.74 

-144*1 . 43 

0.010 

0.040 

0.150 

0.040 

ISC .703 

1.500 

23B.G00 

0.874 

0.SS3 

14613.70 

-7270.37 

o.oio 

0 . 053 

0.137 

0.053 

157.703 

2.000 

317.394 

1.860 

■ C.E98 

11797.40 

-4350.46 

0.010 

0.065 

0.125 

0.065 

156.704 

2.500 

395.402 

3.350 

1.311 

1 0 0 f 9 . 33 

-2098,75 

0 .010 

0.075 

0.115 

0.073 

15S.704 

3.000 

473.055 

5.443 

1 .726 

. 8795.53 

-2U65. 07 

0 .010 

0.035 

0.105 

©.DCS 

1 5 4* . 7 0 5 

3.500 

550.108 

8.201 

2.317 

7902.15 

-1545.03 

0.010 

0 . 073 

0.097 

0.093 

153.705 

4.000 

626.631 

11 .681 

2.B99 

7221.36 

-1199.32 

0.010 

0.100 

0.090 

0.100 

1S2.706 

4.50 0 

7C2.614 

15.940 

3.527 

662S.27 

-950.19 

0 .010 

0.107 

0 . 033 

0.107 

151 .707 

5.000 

773.045 

21.028 

4.193 

62S1 .96 

-783.57 

0.010 

0.113 

0.077 

0,113 

ISC. 707 

5.500 

052. 911 

26.983 

4.907 

5894.20 

-653.27 

0.010 

0.118 

0.072 

0.110 

149.708 

6. 00 0 

927.193 

33.843 

5.652 

5593.50 

-553.62 

0.010 

0 .123 

0.067 

0.123 

146.708 

- .500 

1000. ESQ 

41.637 

6.430 

5336.06 

-47S.84 

O.Olt 

0.127 

0.063 

0.127 

147.70? 

i .00 0 

1073.971 

50.309 

7.237 

SI 15 . 72 

-410.29 

0.011 

0.131 

0.059 

0.131 

146.709 

7.S0C 

1146. 424 

60.122 

8. C72 

4724 .:-6 

-357.43 

0.011 

0.134 

0.055 

0 . 134 

145.710 

8.000 

1218. 231 

70.052 

G .932 

4756.50 

-315.01 

0.012 

0.137 

0.051 

0.137 , 

. 144.710 

8.500 

1209.373 

02.591 

9.814 

4607.90 

-200.50 

0.012 

0.13? 

0.048 

0.13? 

143.711 

9.000 

1359. C32 

9S.34B 

10.718 

4474.95 

-252.12 

0,013 

0.141 

0.046 

0.141 

142.711 

9.500 

1429. S80 

109.130 

11.641 

4354.95 

-220.54 

0.013 

0.143 

0 .044 

0.143 

141 .712 

10 .000 

1493.622 

123.940 

12.592 

4243.75 

-200.7? 

0.014 

0.145 

0 . 042 

0.145 

14C .713 

1C. 501 

1566.914 

139.779 

13.539 

4145.63 

-192.11 

0.014 

0 .146 

0 .040 

0.146 

13°. 713 

11.001 

1634.445 

156.645 

14.512 

4053.21 

-177.93 

0.014 

0.147 

0.038 

0.147 

133.714 

11.501 

1701.195 

174.535 

15.500 

3967.34 

-165.81 

0.014 

0. 14B 

0.037 

0.146 

137.714 

12.001 

1767.144 

193.443 

16.501 

3607.11 

-153.30 

0.015 

0 .14? 

0.036 

0.149 

136.715 

12.501 

1G32.273 

213.361 

17.515 

3C\1 .72 

-146.35 

0.0 V5 

0.150 

0.035 

0.150 

135.715 

13.003 

1896.561 

234.280 

10.5*0 

3740.55 

-130.52 

0.015 

0.151 

0.034 

0.151 

134.716 

13.502 

1959.990 

256.191 

19.570 

3673.03 

-131.67 

0.015 

0.151 

0.033 

0.151 

133.716 

14.003 

2022.541 

279.000 

20.626 

5*08.73 

-123.63 

0 . n 16 

0 . 152 

0.033 

0.152 

132.717 

14.502 

2084.194 

302.935 

21 . 6SS 

3S47.23 

-120.40 

0.016 

0 . 152 

0.032 

0.152 

13: .718 

15.002 

2144.930 

327.742 

22.753 

3403.22 

-115.73 

0.016 

0.152 

0.032 

0.152 

130.710 

15.503 

2204.733 

353. 4B5 

23.832 

3-31 .41 

-111.5? 

0.016 

0.152 

0 . 031 

0.152 

129.71? 

16.003 

2263.583 

380 .147 

24.919 

3376.56 

-107.90 

0 . 0i6 

0.152 

0.031 

0.152 

12S.719 

16.503 

2221.463 

407.711 

26. C J6 

3323.44 

-1 04.60 

0.017 

0.152 

0.031 

0 . 152 

i2'.7?n 

17.002 

2370.355 

436.15* 

27.121 

32 > 1 . 09 

-101.64 

0.017 

0 . 1 52 

0.031 

0 . 15? 

12; .?20 

17.502 

2434.245 

465. 47 U 

23.235 

z. : 1 . 75 

— 9G . 

0.017 

0 . 153 

0.P31 

0.152 

12- . 721 

10.002 

2489.113 

495.626 

29.337 

3172.07 

-96.57 

0.017 

0 . 1S2 

0.031 

0.152 

124.721 

IS. 501 

2542.946 

526.605 

30.487 

3125.14 

-94.30 

0.017 

0.152 

0.031 

0.152 

122.722 

19.001 

2595.725 

55G.305 

31 .625 

3078.46 

-92.33 

0 .oic 

0,152 

0.031 

0.15? 

122.722 

17.501 

2647.439 

590.945 

32.771 

3032.73 

-90.56 

, 0.018 

0.152 

0.031 

0.152 

121.723 

FINAL 1 

VALUES 










19.639 

2661.522 

600.065 

32.771 

3020.26 

-90.09 

0.010 

0.152 

0.031 

0.152 

121 .447 
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TABLE 20 


TERPS CATEGORY I) - CONCORD 

NASA HIGH SPEED EXIT ANALYSIS 
INPUT SUMMARY 

PROBABILISTIC MODEL 


AVERAGE STD DEV UNITS 


TOUCHDOWN LOCATION 
TOUCHDOWN SPEED 
TIME-TOUCHDOWN TO DECEL 
DECELERATION RATE 
50% PROBABILITY SPEED 
DIST: THRESHOLD-START OF TURN 
SPEED AT MAJOR DECEL REDUCE 
NOSE-WINGTIP DIST 


1500 . 

30 . 

FEET 

275. 

15. 

FT/SEC 

5. 

1 . 

SECONDS 

6.00 

1.00 

FT/SEC** 

201 . 

i. 

FT/SEC 

6300 . 


FEET 

150 . 


FT/SEC 

150. 


FEET 


AIRCRAFT TRACKING PROGRAM 


AIRCRAFT CHARACTERISTICS 

WEIGHT (LBS) YAW INERTIA 

245000.0 12664509.0 

WHEELBASE 

59.7 

% LOAD ON MAIN GE^R 
94.7 

FLIGHT CONDITIONAL 
EXIT VELOCITY 
185.009 

PARAMETERS 
MAX EXIT ANGLE 
37.000 

MAX Y MU 
0.200 

DECEL RATE 
-2.000 

SYSTEM WIDTHS 

RUNWAY 
150 .00 

WINGSPAN 

83.80 

CLEAR DIST. 
116.90 

DIST TO TAXI WY 
600.00 

SCRUB CHARACTERISTICS 
Y MU RADIUS 




0.0000 0.1E+32 
0.0100 1000 . 0 
0.0200 000.0 
0.0300 50.0 
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TABLE 21 


CALCULATED RESULTS 


CRITICAL VALUER: 

ascs*-* a/.c MinF roir.r 

iOr >t f;. 168.60 0.57 6„».OO.UU 


xmm 

I'M \ .72 


i p: r.r 
i.v;» 


PROD 

AVG TIME 

STD DEV 

XT I ME 

ttime 

YA 

1.0000 

23.0677 

1 .423? 

31.4201 

39. V709 

201 . 0000** 

0 .9V99 

25.B&77 

1.4237 

31.4201 

39.9909 

2CO.COOO** 


CONSTANT NLG SIDE LOAD ANALYSIS 
PASSENGER COMFORT IGNORED 


TIME 

X 

Y 

THETA 

RADIUS 

R RATE 

MU SCRUB 

MU CENT 

HU I 

C CENT 

VELCCI 

0.000 

0 .000 

0.000 

0 . OOP 

100000 .DO 

-559291 , 56 

0 .010 

0.011 

0.179 

0.011 

185.80? 

o .soo 

92.6S4 

0.001 

0 . 124 

27052.62 

-36965. 2C 

0.010 

0.030 

0.1S2 

0 . 038 

184 .B1C 

1 . o n o 

184 .BOO 

0.464 

0.373 

17249.42 

-17123.59 

0.010 

0.061 

0 .129 

0 .C61 

103 .310 

1 .soo 

276 . 4SB 

1.331 

0.723 

13007.73 

-5915.4- 

0.010 

o.cso 

0.110 

0.080 

182.811 

2.000 

367.601 

2.831 

1 .173 

1C740 .90 

-3473.17 

0,010 

0 . C 76 

0.094 

0.096 

101 .611 

2.56 0 

4 SC . 227 

S.0Q7 

1.693 

9335,67 

-7273.05 

0.010 

0.109 

0.081 

0 .1C? 

100,612 

3.000 

548 . 329 

8.203 

r t 27 ? 

0334.40 

-159 l . Pc 

C.010 

0 . 120 

0.070 

0.120 

179.812 

3. SOO 

637.092 

12.260 

2.91? 

7700.65 

-1 1/3.64 

o . oio 

0.129 

0.061 

0.129 

170.313 

4.000 

726.904 

17.327 

3.6C7 

7186.90 

-053.74 

0.010 

0.137 

0 . 053 

0.137 

177.813 

4. SOO 

BIS. 347 

23.460 

4.325 

6707. SO 

-709.26 

0.010 

0 .143 

0.047 

0.143 

176.014 

S. 000 

903. 2CS 

30.702 

5.057 

6469. S4 

-573.83 

0.010 

0.143 

0.042 

0.146 

175. 81S 

S.SOO 

990.459 

39 . 009 

S.OVO 

6207.70 

-474 23 

0.010 

0 .153 

0.C37 

0.153 

174.015 

a. oon 

1 0 77 . C 89 

48 . 643 

6.710 

SWO .14 

-379.26 

0 . CIO 

0.157 

0.033 

0.157 

173.016 

..soo 

1163.077 

59.400 

7.552 

SODS. 46 

-341.82 

C.010 

0 .160 

0 . 030 

C .160 

172.816 

/. ono 

1248. 399 

71.362 

e . 4i4 

S646. 16 

-297. 1C 

0.010 

0 .162 

0.028 

0.162 

171.817 

7. soo 

1333.036 

□4.543 

9.293 

SS06.75 

-261.34 

0.010 

0 . 165 

0.025 

0.163 

170.617 

8.000 

1416. 96S 

90.950 

10 . 151 

5303 . 1 0 

-233.72 

0.010 

0 .166 

0.024 

0 .166 

165.818 

B.S00 

1S00 . 167 

114.506 

11.102 

5272 . OO 

-211 . 10 

0 . 010 

0.169 

0.022 

0.166 

160.618 

9.000 

1S82.617 

131 . 452 

12.025 

5171 .27 

-192.75 

0.010 

C .169 

0.C21 

0.16? 

167. Gl? 

9. SOO 

1664. 29S 

149.546 

12.961 

5078.76 

-177.74 

0.010 

0.170 

0.020 

0. 17C 

166.819 

10.000 

1745.179 

168. 062 

13.907 

4793 . 07 

-165. 4C 

0.010 

0.171 

0.019 

0.171 

165. e2G 

10 .SOI 

1C25.246 

169.394 

14.063 

4913.00 

-155. IE 

0.010 

0 .172 

0.018 

0 . 172 

164.821 

11.001 

1904.477 

211.134 

15.829 

4037.57 

-1^5.67 

0.010 

0 . 17Z 

0.018 

0 . 172 

163.021 

11 .SOI 

1902.848 . 

234 . 072 

16.CC3 

4766. C? 

- i 3 ~ . 5 4 

0.010 

0.173 

0.017 

0.173 

162.922 

12.001 

2060.339 

250.196 

17.756 

4697 .06 

125. V J 

o . o i t 

0.173 

0.017 

C. 173 

161 . 822 

12. SOI 

2136.929 

233.495 

13.777 

4632.40 

-123.44 

0.010 

0.173 

0.017 

0.173 

160.223 

13.002 

2212. 59S 

309. V54 

19.775 

4569.29 

-124 . 1 1 

m . Jin 

0 . 174 

0.016 

0 .174 

159. 823 

13.S02 

2287.319 

337.559 

23.701 

4508 . 18 

-120.35 

0.010 

0 . 174 

0.016 

0.174 

1 SB , 024 

14.002 

2361 . 080 

366.295 

21.793 

4449. Bl 

-117.16 

0 .010 

0.174 

0.016 

0.174 

1S7.D24 

14 .S02 

2433. B57 

396.145 

22.013 

4390.93 

-114. 4C 

0.010 

0.174 

0.016 

0.174 

156. 82S 

15.002 

2505.631 

427.092 

23.043 

4334.35 

—111. 96 

0 .010 

C . 174 

C.016 

0.17* 

1S5.G26 

IS .503 

2576.332 

459.117 

24 .b; 3 

4272.92 

-tc o .r:' 

0.010 

C . 174 

0.016 

0.174 

154.825 

1 & . 003 

2646. 09C 

492.703 

25.912 

4224 . 49 

-107. ?.? 

0.010 

0.174 

0.016 

0 . 174 

153.027 

16 .S03 

2714.737 

526.330 

25. 

4170.96 

21 

o.oin 

0.174 

0.016 

0 . 174 

152.027 

17.002 

2782.304 

561 .470 

73. C*. 2 

4ii5. :s 

- 1 C 1 . 6*2 

0.0)0 

0.174 

0.016 

0.174 

1S1 .820 

17.502 

21)40.772 

577 . 627 

2V.07; 

4066 . 27 

-103.;.’- 

0 . 01C 

0.174 

0.016 

0 . 174 

ISO. DEE 

FINAL 

VALUES 










17.S3S 

2053.121 

600.047 

29.071 

4C62.06 

-1D3.2C 

0.0)0 

0 .174 

0.016 

0 . 174 

150.762 
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TABLE 22 


TERPS CATEGORY E - LOCKHEED F-104 

NASA NIGH SPEED EXIT ANALYSIS 
INPUT SUMMARY 

PROBABILISTIC MODEL 


AVERAGE STD DEV UNITS 


TOUCHDOWN LOCATION 
TOUCHDOWN SPEED 
TIME-TOUCHDOWN TO DECEL 
DECELERATION RATE 
502 PROBABILITY SPEED 
DIST: THRESHOLD-START OF TURN 
SPEED AT MAJOR DECEL REDUCE 
NOSE-UINGTIP DIST 


1500 . 

30 . 

FEET 

370 . 

25. 

FT/SEC 

5. 

1 . 

SECONDS 

6.00 

2.00 

FT/SEC*#2 

292. 

1 . 

FT/SEC 

8400 , 


FEET 

150 . 


FT/SEC 

35. 


FEET 


AIRCRAFT TRACKING PROGRAM 


AIRCRAFT CHARACTERISTICS 

WEIGHT (LBS) YAW INERTIA 
22000.0 137776.0 


WHEELBASE 
IS. 0 


FLIGHT CONDITIONAL PARAMETERS 


EXIT VELOCITY 
276 . 225 


MAX EXIT ANGLE 

20.000 


MAX Y MU 

0.200 


SYSTEM WIDTHS 

RUNWAY WINGSPAN CLEAR DIST. 

150.00 21.90 85.95 


SCRUB CHARACTERISTICS 
Y MU RADIUS 

0.0000 0.1E+32 

0.0200 50.0 


7 . LOAD ON MAIN GEAR 
92.3 


DECEL RATE 

- 2.000 


DIST TO TAXIWAY 
600 . 00 
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TABLE 23 


. mLCULAU P KLUJLIfi 
CRITICAL VALUES: 


A5CSV 

276.22 

ASC 

. 264.01 

PTINL 

6.06 

CP 1ST 
0400. 00 

XD1ST 

1634.71 

I /I SI 
10034.71 

PROD 

AVG TIhC STD DEV 

XTim 

T T 1 ML 

T A 


1.0000 
0 ;99?9 

2-1.7617 2.3551 

24.7617 2.3551 

33.9466 

33.9466 

40.0104 

40.0104 

292. 00004* 
271.00004* 



CONSTANT NL5 SI DC LOAD ANALYSIS 
P ASSENCcR COM OnT IGNORED 


TIME 

X 

Y 

THETA 

RADIUS 

R RATE 

MU SCRUB 

MU CENT 

MU I 

C CENT 

VELOCITY 

0.000 

0.000 

0.000 

0 .000 

lonono . co 

-963755.37 

0 . 

0 . 024 

0 . 156 

0 . 024 

276 . 225 

0 .soo 

137.060 

0.223 

0.231 

22530 .26 

-23733 . 96 

0.020 

0.104 

0 . 076 

0.1C4 

275.217 

1.000 

275 .212 

1 .256 

0.654 

16352.71 

“6175.02 

0 . 020 

0 .143 

0.037 

0.143 

274 .210 

1 .50 0 

412.047 

3.419 

1 .170 

14333.35 

“2428.62 

0 . 020 

0 . 161 

0.019 

0.161 

273 .203 

2.000 

540 . 353 

6,062 

1 .731 

13537.84 

-1150.20- 

0 .020 

0.170 

0.010 

0.170 

272 . 176 

2. SOO 

68 4 , 1 1 4 

11.655 

2.316 

11111.43 

“626.37 

0 . 020 

0 . 174 

0.006 

0.174 

271 . 1 P? 
270 . 152 
269. 175 

3.000 

019.314 

17.827 

2.914 

12864 . 03 

“373 . 05 

0.020 

0.176 

0.004 

0 . 176 

3. SOO 

VS3.V40 

25.309 

3.519 

12697. e9 

-204 . 74 

0.020 

0 . 177 

0 .003 

0,177 

4.000 

1087.976 

34.344 

4,128 

12569.92 

-233.25 

0.020 

0 .178 

0 .002 

0.178 

260. 168 
267.161 
266. 154 
265 . 147 

4. SOO 

1221 .407 

44,60? 

4.741 

12460 .26 

-208.30 

0 . 020 

0 . 178 

0.092 

0 .178 
0.173 

S.000 

1354 .217 

56.419 

5.356 

12359.55 

-175.94 

0 . 020 

0 . 173 

0.002 

s.soo 

1406.393 

69.529 

5.974 

12263.33 

-107.63 

0 . 020 

0 . 173 

V .002 

0 . 17S 

6.000 

1617.919 

84.012 

6.595 

12169.46 

-106.20 

0 . 020 

0 . 173 

9.002 

0.178 

264 . 14C 
263 .123 

6. SCO 

1748.779 

99.061 

7.210 

12076.90 

-104.17 

0 .020 

0 ,173 

0 . 0C2 

0.173 

7.000 

1070.960 

1 17.063 

7.043 

11905.21 

-102.07 

0.020 

0 . 178 

0.002 

0.173 

262 . 1 2S 

7.500 

2008.446 

135.627 

8.471 

11374 . C7 

-101 .87 

0 . 020 

0 . 173 

9 .002 

0 . 176 

261 .118 
260.111 
237 . 104 
2SB . C97 
257 . 090 

S. 000 

2137.224 

155.520 

7. 101 

110C3.30 

-131.01 

0 . 020 

0 .178 

9.002 

0 . 17C 

3.500 

2265.277 

176 .76S 

9.734 

11713.05 

-1O0 . 17 

0 . 020 

0 .178 

0.002 

0.170 

9 . CCO 

2372.573 

199.329 

10.369 

11623.21 

-179.52 

0.070 

0. 17B 

0 .002 

0 . 1*76 

7. SOO 

2519 . 155 

223.212 

11 .006 

11533.70 

-178.03 

0 .020 

0.173 

0.002 

0.178 

10.000 

2644.949 

240. 4CiS 

11 .646 

11444 . 49 

-178.09 

0.020 

0 .178 

0 .002 

0 .175 

236 . ZUZ 

10. SOI 

2769.962 

274.900 

12.200 

11355. B2 

-176.77 

0.020 

0 . 178 

0 .002 

0.178 

255 .003 

11.001 

2Cy4 .103 

302. 6GG 

12.933 

1 1267.75 

-175.71 

0.020 

0.178 

0.002 

0.173 

254 .023 

11 .SOI 

3017.576 

331 .760 

13.500 

1 11G0 . 16 

-174.86 

0 . 020 

0 . 178 

0 .002 

0 . 178 

253 . 024 

12.001 

3140 .196 

362.107 

14.230 

11092.93 

-174.08 

0 . 020 

0 .178 

0.002 

0.170 

232 . 065 

12. SOI 

3261 .940 

393.719 

14.002 

11006. 12 

-173.29 

0 .020 

0 .170 

0.002 

0.173 

251 . 0G5 
250 . 0 06 

13.002 

37.02.843 

426. 5B7 

15.537 

10919.70 

-172.61 

0 . 020 

0.170 

0.002 

0.178 
0 .178 

13.502 

3502.082 

460 .701 

16.194 

10333.64 

-171.95 

0 . 020 

0 .178 

0.002 

24?, 036 

14 . C 02 

3622.043 1 

476.050 

16.854 

10747.07 

-171 . 23 

0.020 

0.178 

0.002 

0.17Q 

240.087 
247 . 037 

14.S02 

3740.307 

532.624 

17.517 

10662.42 

-170.47 

0.020 

0 . 178 

0.002 

o . i ; c 

15.002 

3857. 667 

S70.413 

10.102 

10577.43 

-169.87 

0.020 

0.178 

1.002 

0.178 

246.088 

FINAL 
IS. 384 

VALUES 

3946.475 

600.017 

10.182 

10512.82 

-167.31 

0.020 

0. 17B 

0.002 

0.170 

245.326 
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Optimum High Speed Exit Paths 


The 747 was used as an example aircraft because it is the largest, and 
therefore most difficult commercial aircraft to maneuver* Designing an exit 
for the 747 would insure that it could be used at the design speed by smaller 
aircraft* 

Figure 7 shows three optimum exit turn paths for the Boeing 747. These paths 
were based on a limiting nose gear side y of 0.2 and an aft c.g. Two speeds 
were used — the 27 m/s (60 MPH) design speed of the present FAA exit and a 96 
knot (110 MPH or 49 m/s) exit speed resulting from the 747 probabilistic 
studies. Additionally, the 27 m/s (60 MPH) exit was determined for both a 
constant-speed exit and a favored 0.61 m/s 2 (2 fpss) deceleration while in 
the exit. The deceleration resulted in reducing the longitudinal length of 
the exit by 16%, from 433 m (1422 ft.) to 364 m (1195 ft.). 


C N - 0.2 LIMITATION 

MAX 
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FIGURE 7. OPTIMUM HIGH-SPEED EXIT PATHS - 747 
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PROGRAM LISTING 

PROBABILISTIC MODEL OF OPTIMAL HI^H SPE P D 
RUNWAY TURNOFFS 

REAL NTIP 

REAL *4 NTlTl.E < 20 ) 

DIMENSION ARDSd6) .DPXd6> .P( 16) .DP (16) 
DIMENSION SND<400) 

ARBT (15) .ARBI (15) . DT ( 15) 

PEX ( 16 ) .ZROTCIS) .2(16) .ZEXd6> 
CALL ASSIGNd . 'HISPEX.DAT' ') 

CALL ASSIGN < 2 \ » SNDtDAT j !»->•-• 

CALL ASSIGN ( 3 . ' SOLVIT . DAT s 1 * ■) 

READ (1 .11) (NTITLE(I) .1=1 .20 ) 

FORMAT < 20A4 ) 


DIMENSION 


10 


READ(1 .10) l)A, 
READ d . 1 0 ) XA, 
F0RMAK6F5 . 0 ) 
IBIRD=0 
IFLAG=0 


US. 

xs'. 


DA. 

VS. 


US 


WA. 

. ZA, 


NTIP 


1.25 


13? 


1 

P 


WRITE(S . 125) 
FORMAT ( ' ENTER 

> 

READ ( 5 > *) IURT 
WRITE (5.13?) 
FORMAT ( * ENTER! 


READ ( 5 > * ) IPRT 
UR I TE ( 5 . 6 ) 
FORMAT ( * ENTER 


5 

6 


FOR 

FOR 


TERMINAL 

HARDCOPY 


PRINTOUT V . 
PRINTOUT ' ) ' 


FOR 

FOR 

FOR 


STANDARD 


PRINT ' ./ 
INTERMEDIATE VALUES' 
FINAL VALUES ONLY ' ) 


i 

p 


FOR 40 SEC OCCUPANCY'/. 

FOR MAXIMUM PROBABILITY'/. 
ENTRY OF YA VALUE' ) 


ENTER INITIAL YA SPEED') 


,1 = 1 , 20 ) 


750 

20 


C 

C 

C 


1 TO SOLVE 
’ 2 TO SOLVE 

' 3 FOR USER 

READ ( 5 . *) IYA 
WRITE (5. 7) 

FORMAT ( ' 

READ(5,*)YA 

IFdPRT . EO . 3) GO TO 750 
WRITE( IURT .3) (NTITLE(I) 

FORMAT (29X .20A4/) 

WRITE ( IWRT , 1 ) 

FORMAT ( 30X , ' NASA HIGH SPEED EXIT ANALYSIS'/ 

135X. ' INPUT SUMMARY'// .T35 . 'PROBABILISTIC MODEL'./ 

2/T43. 'AVERAGE STD DEV UNITS'/) 

WRITE ( IWRT ,2>UA .US .VA, VS .WA .US .XA . XS . YA . YS . C . ZA . NTIP 
FORMAT d OX , 'TOUCHDOWN LOCATION' .T40 .2F10 . 0 ‘ T65 . ' FEET ' 
i/i OX .' TOUCHDOWN SPEED '. T4 0 . 2F1 0 . 0 . T65 .' FT/SEC '. /j OX 'TIME-' 
2. 'TOUCHDOWN TO DECEL' ,T4 0 .2F10 . 0 . T65 . 'SECONDS' ,/i0X 
2 'DECELERATION RATE ' .. T40 ,2Fi 0 . 2 , T 65 , ' FT/SEC**2 ' . 

3/10X.. '50% PROBABILITY SPEED ', T4 0 , 2F 1 0 . 0 . T65 , ' FT/SEC ' /IflY 
3 ' DIST : THRESHOLD-START OF TURN ' . T40 . FI 0 . 6 . T65 . ' FEET ' , 

4/10X . 'SPEED AT MAJOR DECEl. REDUCE'. T40.F10.0.T65. 'FT/SEC'/ 

510X. 'NOSE-UINGTIP DIST' .T40 .FID . 0 .T65. 'FEET'//) ‘ 

READ (2.20) SND 
FORMAT ( 1 0F5 . 4 ) 

YAPERM=YA 


CALCULATE DISTANCES 
A = UA + VA*WA 

B = A ^ (VAs*2-ZA**2)/(2.*XA) 
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c 

c 

c 


021 

825 

600 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

500 


C 

C 

C 


C 

C 

C 


100 

110 


CALCULATE SPEEDS DURING LANDING 

ASH = (VA*#2-2*XA*(B-A) )#*0 , 5 
IF(VA**2.GT. (2*XA*<C-A) ) ) GO TO 825 
URITE< IWRT 021 ) 

FORMAT < ' SPEED AT EXIT IS NEGATIVE') 

STOP 

CONTINUE 

ICK=i 

FORMAT ( ' AT POSITION '.12) 

IF<IPRT.E0.2)WRITE(5.600)ICK 

ICK=2 

ASC = <VA#*2-2#XA*<C-A> )**0 .5 
ASCSV = ASC 

CALL EXPATH TO DETERMINE EXIT PATH 
TIME AND DISTANCE 

CALL EXPATH ( ASC , PTIME.XDIST. IURT_. IPRT. I FLAG ) 

UPON RETURNING FROM EXPATH. THE VALUE OF ASC 
<AVG SPEED AT START OF TURN) HAS DEEN REDEFINED 
TO SPEED AT END OF TURN. AFTER THIS POINT. 

AVG SPEED AT START OF TURN HAS BEEN REDEFINED AS 'ASCSV' 

TDIST = C + XDIST 

CONTINUE 

MS=' » 

MSP=' ' 

ASC2=YAPERM- < 1 . S*YS ) 

IF (ASCSV . LT .ASC2)MS='#' 

IF < IPRT . EO .2)WRITE<5.600) ICK 
ICK=3 


CALCULATE STANDARD DEVIATIONS 


SDA = <VS**2+(US*XA/VA>**2 + 
IF ( IPRT . EO .2)URITE(S. 600)1 CK 
ICK=4 


<US*XA)**2)**0 .5 


SDB = <SDA**2 + <2*XS*(B-A)/(VA+ZA>)**?)**0 

IF(IPRT . EO . 2 ) WR ITE (5 . 60 0 ) ICK 

ICK=5 


.5 


SDC - <SDA**2 + < 2*XS* ( C-A ) / < VA+ASCSV) ) **2) **0 

IF( IPRT. EO. 2) WRITE (5 .600) ICK 

ICK=6 


5 


CALCULATE OCCUPANCY TIMES 
ASA = VA 

TA = ( (UA+NTIP)/VA) + WA 

IF (B.GT.C) GO TO 100 

TB = TA + 2*<B-A)/<ASA+ASB) 

TC = TB + 2* < C-B) / < ASB+ASCSV ) 
GO TO 110 

TC = TA 4 2#(C-A)/< ASA+ ASCSV ) 
CONTINUE 
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CALCULATE DATA FOR AIRCRAFT TRAVELING 
ONE STANDARD DEVIATION BELOW NORM 

CALCULATE DISTANCES 

AP = A 

DP = AP + <D-A)>k( < ASA-SDA ) -ZA > / ( ( ASA— SDA ) - ( ZA— SDB) ) 
CP = c 

CALCULATE SPEEDS 

ASAP = ASA - SDA 
' ASDP = ZA 

ASCP = ASCSV - SDC 
ASCPi = YA - 1.5*YS 
IF( ASCP i .GT. ASCP)MSP='*' 

CALCULATE OCCUPANCY TIMES 

TAP = 2*A/(VA-VS+ASAP) 

IF (BP.GT.CP) GO TO 120 

TDP = TAP + 2*<BP-AP)/!ASAP+ASPP> 

TCP = TDP + 2*(CP-BP>/!ASBP+ASCP> 

GO TO 130 

120 TCP = TAP + 2*!CP~AP ) / < ASAP+ASCP ) 

130 SDOT = TCP - TC 

IF ! SDOT . GT . 0 . ) GO TO 900 
WRITE (5.901) 

901 FORMAT ( ' STANDARD DEVIATION LESS THAN ZERO') 

STOP 

900 CONTINUE 


WRITE INTERMEDIATE VALUES FOR CHECKOUT 


135 


136 


IF( IPRT . NE . 2 ) GO TO 138 

WR ITE ! I WRT , 1 35 ) A , D , ASB . ASC , SDA . SDB , SDC . ASA . TA . TB TC 
FORMAT ( ' INTERMEDIATE RESULTS IN ORDER:'/, ’ 

^ * A , D , ASD , ASC , SDA , SDE< , SDC . ASA , TA ,TB . TC ' / . ' AP . BP CP ASAP ASBP^ 
2' ASCP .TAP TAP , TBP , TCP , SDOT ' , // 1 1F8 . 2) ' ' ' ' * 

WRITE! I WRT , 136) AP , BP .CP , ASAP , ASBP .ASCP .TAP .TBP .TCP .SDOT 
FORMAT (10F8. 2) ' ... ,ouui 


CALCULATE PERCENT AIRCRAFT EXITING 

SPECIFY ARBITRARY SPEED RANGE 

138 CONTINUE 

ITEMP = ASCSV-3 . 

ARBS(i) = ITEMP 
DO 140 1=2,16 

140 ARBS! I ) = ARBS ( I- 1 ) + 1. 

C 

C CALCULATE PROBABILITY FOR SPEED RANGE 
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ISO 

c 

c 


ZZ=(ANS<ZU> > + C, 01) *100. 

ICNT=ZZ 

P d ) = 0 .5+SNDdCNT) 

IF(ZCI>.LT.O. )P<I)=0. 5-SND ( ICNT ) 
CONTINUE 


CALCULATE- DISCRETE 


INTERVAL 



DP ( 1 ) = P(i) 

TOT = P(i) 

DO 160 1=2.15 

dp ( i ) = p (i> - p<i-D 

160 TOT = TOT + DP ( I ) 

DP ( 16) = 1 . -TOT 
C 

C CALCULATE PROBABILITY OF EXITING 

C 

DO 170 1=1,16 

ZEX(I) = <YA-(ARBS< I )-0 , 5) )/YS 
IFIZEXd) .GT.3.99) ZEX(I>=3.9? 
IF(ZEXd) .LT.-3.99) ZEXd>=-3.99 
ZZ=(ABS(ZEXd)> +0.01 >*100 
rcNT=zz 

F’EX (I ) = 0 . 5+SNDdCNT > 

IF(ZEX(I) .LT.0. )PEX(I)=0.S-SND<ICNT) 
170 CONTINUE 
C 

C CALCULATE DISCRETE PROBABILITY AND SUM 

W 

PX = 0. 

DO '180 1 = 1.16 
DPXd) = DP d )*PEX( I ) 

180 • PX = PX + DPX(I) 

c 

C WRITE INTERMEDIATE VALUES 

C 


C 

C 

c 


IF(IPRT.NE.2) GO TO 185 
WRITEdURT. 181 ) 

l81 « T f'ORMAT(/Ti0d PERCENT AIRCRAFT EXITING CALC'/ > ARFtq , 

WRITE < IURT. 183) DP (I ) .ZEX(I) .PEX(I) .DPX(I) 

UR ITE (IURT . 107 > ARBS < I ) . Z d ) .P(I) ‘ 

FORMAT ( 24X . F8 . 4 , F8 . 3 . 2FB .A) 

FORMAT (F8.1 .F0.3.F8.4) 

UR ITE dURT . 188 ) PX 

FORMAT (/Ti 6 ' PERCENT AIRCRAFT EXITING = ' T43 F10 41 
CONTINUE ,.i i («l.H0.4) 

CALCULATE MINIMUM OCCUPANCY TIME 


184 
183 

187 

188 

185 


145 

147 


PXC=ABS ( PX-0 . 5 > 

DO 145 1=1.400 
IFdPXC.GT.SNDd) ) 
URITE<5. 147) 


.AND. (PXC.LE. SND( 1+1))) GO TO 146 


FORMAT (' PX VALUE NOT FOUND') 
STOP 
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146 


BAN=I 

PXCH=BAN/100 . 

0TM1N = TC - PXCH * SDOT 
IF ( IPRT . EQ.2)URITE(IWRT.19i)0TMIN 

191 FORMAT CT10 . * MINIMUN OCCUPANCY TIME= » , T43 , Fi 0 . 2 . » SEC») 
C 

C CALCULATE AVERAGE RUNWAY OCCUPANCY TIHE 

C •••• •• 

C 

C INITIALIZE ARBITRARY OCCUPANCY TIMES 

C 

ITMIN=OTMIN 

ARBT(i) = ITMIN 

DO 190 1=2. IS 

ARBT(I) = ARBT(I-i) + 1. 

190 CONTINUE 

ARBT ( 1 > =OTMIN 

COMPUTE Z VALUES 

DO 200 1 = 1. IS 

ZROT(I) = (ARBT< I )-TC)/SDGT 
IF(ABS(ZROT(I>) .GT.3.99) GO TO 201 
ZZ=(ABS(ZROT(I))+0.0i)*100. 

ICNT=ZZ 

P < I > = 0 . 5+SND ( ICNT) 

IF(ZROTd) .LT.O . >P ( I >=0 . 5-SND( ICN V > 

200 CONTINUE 

201 IBAN=I-i 

CALCULATE INTERVAL MIDPOINTS 
ARBI ( 1 ) = (ARBT(l)+ARBT(2)>/2. 

• DO 210 1=2. IRAN 
210 ARBI ( I ) = ARBT ( I ) + 0.S 

CALCULATE DISCRETE PROBABILITY AND SUM 

PSUM=0 . 

TSUM=0 . 

P ( IDAN+1 ) =1 . 

DO 220 1=1 .IBAN 

dp ( i ) = p(i+i> - p(i) 

DT ( I ) = DP(I)*ARBI(I> 

PSUM = DP ( I ) + PSUM 
TSUH = DT ( I ) + TSUM 
220 CONTINUE 

AROTEX = TSUM /PSUM 

WRITE FINAL PARAMETERS 

IFUPRT.NE.2) GO TO 230 
WRITE ( IWRT .231 ) 

231 FORMAT ( / > AVERAGE RUNWAY OCCUPANCY TIME CALC'/ 

1 * AR BT » . T 1 2 . » ZR OT ' , T2 1 , ' P * . T29 , >m B I ' . T37 . ' DP ' ! T4S . ' DT ' / ) 

DO 232 1=1. IBAN ' ' ’ * 

WRITEt IWRT .233) ARBT ( I > .ZROTCI) ,P(I) 

232 UR ITE ( IWRT .234 ) ARBI ( I ) . DP ( I ) . DT ( I ) 
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233 FORMAT <F0.2.F8.3,F8.4> 

234 FORMAT ( 24X . F8 . 2 . F0 . 4 , FO . 3 ) 

WRITEdURT, 235>PSUM,TSUM 

235 FORMAT (/' AREA SUM= * . T20 . Fi 0 . 4 . / ' TIME SUM= ' . T20 , FI 0 . 2 > 

. .... WRITEdURT, 238>AR0TEX ' 

23B^'^vFOJlH^Oifl > v..!-i:AyERAGE RUNWAY OCCUPANCY TIME=» .T43.F10 SEr»i 

230 CONTINUE ' - ~ X ~ 

CALCULATE YTIME ANDTTIME 

XTIME=AR0TEX+3 . ?)KSI)OT 
TTIME = XT I ME + PTIME 
IFdPRT . EO . 2 > IJR ITE ( I WRT , 65 i ) SDOT 
651 FORMAT di 0 , ' STANDARD DEVIATION = ' . T43 , FI 0 . 2 . > SEC') 

IF < IPRT . EO . 2) WRITE d WRT , 650 )XTIME ' 

650 FORMAT di 0 . ' AVG TIME PLUS 3 STD DEV= ' . T43 . Fi 0 2 * SEC') 
URITE<3,700) PX, AROTEX. SDOT. XTIME, TTIME. YA. MS. MSP 
IF ( IYA . EO . 3 > WRITE <5. .1327) PX ' • • • 

827 FORMAT ( ' PROS = ',F7.4) 

IF< (XTIME .LT .39 .99) . OR . (XTIME .GT . 40 . 01 ) ) GO TO 800 
810 CONTINUE 

PXX=PX*i00. 

. WRITE( IWRT .853)PXX 

853 FORMAT(/T10,» PERCENT AIRCRAFT EXITING = '.T43.F10.4) 

WRITEdWRT ,238)AR0TEX 
WRITEdWRT, 651) SDOT 
WRITEdWRT, 650) XTIME 
WRITE (I WRT. .850) YA 

850 FORMAT (T10 , ' FINAL 50% PROBABILITY SPEED= ' . T43 .FI 0 . 2 . > M/SEC') 

WRITEdWRT. 852) ASCSV ... 

852 FORMAT (T1 0 , ' AVERAGE SPEED AT EXIT=' ,T43 .F10 . 2 . » M/SEC' > 

WRITEdWRT. 851) IBIRD 

851 FORMAT (/T20,' PERFORMED IN '.13.' ITERATIONS V) 

GO TO 801 

800 CONTINUE 

IF(IYA.EO.l) YA=YA*XTIME/40 . 

IF (IYA . EO ,2) YA=YA-i . 

IF(IYA.LT.3) GO TO 824 
WRITE <5. 823) 

823 FORMAT (' ENTER NEW YA SPEED: 999 TO END * ) 

READ(5.*)YA 

IF ( YA . EG . 999 . ) GO TO 801 
700 FORMAT ( 6F 12.4. 2Ai ) 

824 IBIRD = IBIRD +1 
IFdBIRD.LT. 15) GO TO 500 

801 CONTINUE 
WRITEdURT. 820) 

820 FORMAT (/' CALCULATED RESULTS'/) 

WRITEdURT, 822) ASCSV, ASC , PTIME , C, XDIST , TDIST 
822 FORMATt/' CRITICAL VALUES :'/ ,T1 0 , 'ASCSV' ,T25 . 'ASC ' .T40 

1 'PTIME ',T55,'CDIST',T70, 'XDIST' . T05 '. ' TDIST ' / . TS , 6 ( Fi 0 . 2 . SX > ) ' 
WRITEdURT, 807) ’ ' 

007 FORMAT ( //6X, 'PROB' , T18, ' AVG TIME' .T30 . 'STD DEV '. T42 .' XTIME ' T54 
1 'TTIME' ,T66, 'YA'/) . ‘ ‘ 

CL05E(UNIT=3) 

CALL ASSIGN ( 3 , ' SOLVIT . DAT : 1 ' ) 

803 RE AD <3. 700 .END=B08)Oi . 02 . 03 . 04 . 05 . 06 . 07 . 00 

UR ITE (I WRT. 700 >01 .02.03.04.05,06.07.08 ‘ 
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GO TO 803 

808 UR ITE ( 5 , 920 ) ' ’ 

920 FORMAT ( / * DO YOU WANT TO CONTINUE EXIT PATH TO TAXI-WAY */ 

1 * 1=YES . 2=N0 . ») m • 

READ ( 5 . * ) I NASA 
IF (INASA .EO. 2) STOP 
- - IFLAG=1 - 

CALL EXPATH (ASCSV . PTIME .XDIST . IURT . IPRT . I FLAG) 

STOP .... 

END 


Aircraft Tracking 
(Subroutine EXPATH) 


SUBROUTINE EXP ATH ( V , PTIME . XDIST . IURT . IPRT . IFLAG ) 
DIMENSION RSCR(20> . YMUSCR(20> 

DIMENSION XP(iOO) .YP(iOO) 

INTEGERS IDO 
REAL JMAX 

CALL ASSIGN (4 . ’PATH . DAT * > 


READ INPUT PARAMETERS 


■ READ(4,10) UT.YAUI .UHLD.PCTM 
READ (4,10) THETMX .YMUMAX . DECEL 
READ (4.10) RUIDTH . USP . TAXDIS 
TUTH = ( RWIDTH+USP ) /2 . 

10 FORMAT (4F12.3) 

DO 20 1=1,20 

20 READ (4,11 ,END=30 ) YMUSCR ( I ) , RSCR ( I ) 

11 FORMAT (F5.4.G10.1) 

30 ICNT=I-i 

REWIND 4 
CLOSE(UNIT=4) 

IF ( IPRT . EQ . 3 . OR . IFLAG . EO . 1 ) GO TO 37 
WRITE( IURT.3S) 


35 


50 


FORMAT ( /T32 , 'AIRCRAFT TRACKING PROGRAM*//) 

UR ITE ( IURT .50 ) UT . YAIJI . UHLD . PCTM 
WRITE ( IURT , A0 ) V, THETMX. YMUMAX. DECEL 
UR ITE ( IURT .75) RUIDTH .USP . TWTH . TAXDIS 
UR ITE ( I WRT ,70 ) 

WRITE (IIJRT, 80) YMUSCR ( 1 ) .RSCR(i) 

URITE(IlJRT,Oi) ( YMUSCR (I).RSCR(I) , I=2.ICNT) 

FORMA! ( T5 , * AIRCRAFT CHARACTERISTICS*/ .Ti 0 . * WEIGHT (LBS ) > 
1 , * YAU INERTIA*, T43, * WHEELBASE* .T57,*X LOAD' ON MAIN GEAR* 
2T10 . 4 ( F12 . 1 .3X)/) 


. T25 

,/ 
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60 FORMAT ( ' * FLIGHT CONDITIONAL PARAMETERS'/ ,T10 . 

1 ' EXIT VELOCITY ' .T2S . ' MAX EXIT ANGLE ' .T45 . ' MAX Y MU'. 

2T 60 . ' DECEL RATE'./ ‘ 

2T 1 0 , 4 ( FI 0 . 3 . 5X ) / ) 

70 FORMAT!/' SCRUB CHARACTERISTICS'/ .TiO . 'Y MU' ,T20 . »RADIU°» /) 

7S FORMAT!' SYSTEM WIDTHS'/ .TI5. 'RUNWAY' .T28 . 'WINGSPAN' . 

tTA2 , 'CLEAR DIST. ' ,T6i , 'DIST TO TAXIWAY VTlo'. 3 (F 10 .2.SXKT61 .FiO 2 ) 

80 FORMAT (TS»F10.4.G10.i) * ‘ ’ ' ' L ' 

81 FORMAT ! T5 .Fi 0 . 4 . F10 . 1 > 

37 JMAX=10 . 0E20 

GMAX=i 0 . 0E20 
INITIALIZE VARIABLES 

85 CONTINUE 

IF ! IF'LAG . EO . i ) THETMX=i80 . 

IF ! IFLAG . EQ . 1 ) TWTH=TAXI)IS 

G=0 . 000001 

T = 0 . 0 

X=0.0 

Y=0 . 0 

THETA=0.0 

R=i 0 0 0 0 0 . 

IPR INT=1 0000 

THETMX=THETMX*3. 1416/180 . 

DT=0 .001 

CONSTANT NLG SIDE MU 


800 

1 


810 


90 

92 


ITR=2 

IF ! I.PRT , EO . 3 . OR . IFLAG . EO . 1 ) GO TO 810 
WRITE (5.800) 

FORMAT!' ENTERS 1 TO PRINT TRACKING RESULTS TO RUNWAY EDGE 
2 TO WAIT FOR PRINT TO TAXIWAY EDGE') 

READ ( 5 . # ) ITR 
CONTINUE 

IF (IFLAG. EQ.l) ITR=i 
IF(ITR.EO.i) WRITE(IWRT.80i0> 

IF < ITR .EO . 1 ) WRITE(IWRT.96> 

IF(ITR.EO.l) WRITE! IWRT . 80 11 ) 

DO 200 IDO=l. 60000 
DO 90 J=2 . 20 * 

IF ( R . GE . RSCR ( J) ) GO TO 92 


t 


/ 


YMUSC=YMUSCR(JJ)-i (YMUSCR! J)-YMUSCR( J J ) ) * ( RSCR ( J J ) -R ) 
/(RSCR(JJ)-RSCR(J)) 

IF ( RSCR (J) . LE . 0 . ) WRITE! IWRT. 8001 ) 

IF(R . GT . RSCR ( 1 ) ) WRITE ( IWRT , 80 0 1 ) 

IF ! RSCR ( ,T ) , LE . 0 .OR. R . GT . RSCR ( 1 ) ) STOP 
YMUC = V**2/(32.2*R> 

YMUI = YMUMAX-YMIJSC-YMUC 

RD = - (YMUI*R**2/(YAWI*V))*WT*WHLB*(PCTM/i00.)*(i.-PCTM/i00 ) 

G = V**2/!R*32.2> ruin/iuu./ 

GDMX = JMAX*(i .-G/GMAX) 

RDMX = —32 . 2*R**2*-GDMX/V**2 
IF(RDMX.GT.RD) RD=RDMX 
IF(Y.GE.TUTH) GO TO 500 
IFdPRINT.LT.500) GO TO 110 
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100 TH = THETA* 180 . /3 . 1416 
IBC=IBC+1 - 
XP(IBC)=X 
YP ( I BC ) =Y 

IF < ITR . EG . I ) UR ITE< IWRT ,8000 ) T.X.Y.TH.R .RD , YMUSC * YMUC . YMUI ,G,V 
IF(RD.EQ.RDMX) WRITE ( IURT . 80 02 ) 

IPRINT = 0 
110 CONTINUE 
C 

C INTEGRATE 

C 

T=T+DT 

V « V+DECEL*DT 
IPRINT = IPRINT + 1 

X = X+V*DT*COS<THETA> 

Y = Y + V£DT#SIN< THETA ) 

THETAD = V/R 

THETA = THETA + THETAD * DT 
R = R+RDfcDT 
200 CONTINUE 
C 

C PRINT FINAL TIME VALUE 

C • 

500 IF (THETA . LE . THETMX) GO TO 520 
WRITE (5*9000) 

9000 FORMAT ( // * WARNING!! WINGTIP NOT CRITICAL POINT! V/) 

STOP 

520 I F < ITR . EG . i > WRITE( IWRT *860 0 ) 

0600 FORMAT <//* FINAL VALUES*/) 

IF ( ITR .EG. i) WRITE<IWRT,8Q00)T,X*Y*TH.R . RD . YMUSC * YMUC . YMUI ,G*V 
PTIME = T 
XDIST=X 
C 

C FORMAT LIST 

C 

96 FORMAT ( */ PASSENGER COMFORT IGNORED * ) 

8000 FORMAT(4F10 . 3 *2F12 . 2 , 5F1 0 .3) 

8001 FORMAT ( 1H *****OUT OF TABLE BOUNDS**** > 

8002 FORMAT ( * LIMITED BY COMFORT* ) 

8010 FORMAT (// * CONSTANT NLG SIDE LOAD ANALYSIS*) 

8011 FORMAT ( //T5 * * TIME * %T17 % * X * ,T27. *Y* ,T35* * THETA* .T46 , * RADIUS * ,T58 
i * R RATE* *T60 * 'MU SCRUB* *T79 * *MU CENT * * T90 * * MU I**T100**G CENT* 

1 .Till, 'VELOCITY*/) 

RETURN 

END 
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